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SUMMARY 


The thermodynamic and transport properties of high-temperature air 
are found in closed form starting from approximate partition functions 
for the major components in air and neglecting all minor components. 

The compressibility, energy, entropy, the specific heats, the speed of 
sound, the coefficients of viscosity and of thermal conductivity, and 
the Prandtl numbers for air are tabulated from 500° to 15,000° K over a 
range of pressure from 0.0001 to 100 atmospheres. The enthalpy of air 
and the mol fractions of the major components of air can easily be found 
from the tabulated values for compressibility and energy. It is predicted 
that the Prandtl number for fully ionized air will become small compared 
to unity, the order of 0.01, and this implies that boundary layers in 
such flow will be very transparent to heat flux. 


INTRODUCTION 


It is axiomatic that the science of aerodynamics must be based on 
a good understanding of the atmospheric medium through which vehicles 
are to fly. Under subsonic flight conditions, air may be treated as an 
ideal gas composed of rigid, rotating diatomic molecules. The thermody- 
namic properties of such a gas are well known and they are accounted for 
in the gas flow equations by the familiar ratio of specific heats, which 
in this case is a constant. Under supersonic flight conditions, air may 
be raised to temperatures where the molecules can no longer be treated 
as simple, rigid rotators. At relatively low supersonic speeds, vibra- 
tional energy is excited and then the specific heats become functions of 
temperature. However, both the thermodynamic and transport properties 
of air in vibrational excitation can be predicted with fair accuracy by 
the methods of quantum statistics and kinetic theory (ref. 1), and the 
air-flow relations can be modified accordingly. Eggers (ref. 2) has 
calculated the effects of vibrational energy on the one-dimensional, 
inviscid flow of diatomic gases, for example. 

Further changes in air properties may occur at still higher flight 
velocity. Flight velocities of practical interest have now increased 
from low supersonic speeds to the escape velocity, 37*000 feet per second. 
Vehicles which travel at these hypervelocities excite the air to such 
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high temperatures that the molecules not only vibrate but may dissociate 
into atoms and even ionize. Under these conditions, the behavior of air 
deviates widely from that of an ideal gas and the thermodynamic and trans 
port properties all become functions of pressure as well as of tempera- 
ture. It is, of course, essential to evaluate these functions in order 
to calculate the pattern of air flow about high-speed vehicles, the vis- 
cous and pressure forces which result, and the heat flux which occurs 
between the air and the vehicle. 

The equilibrium thermodynamic properties of a gas can be calculated 
with great confidence, provided the energy levels of the gas particles 
and the degeneracy of these levels are known. For monatomic and diatomic 
gases this information can generally be deduced from spectroscopic data 
with such accuracy that the calculated thermodynamic properties can be 
trusted to very high temperatures, even where experimental confirmation 
is lacking. In the case of air, however, one of the important energy 
terms was not known with confidence until recently, namely, the dissocia- 
tion energy of molecular nitrogen. This uncertainty arose because the 
available spectroscopic data were consistent with two different models 
for nitrogen dissociation, one leading to a dissociation energy of 7*37 
electron volts per molecule and the other to 9*7 6 electron volts per 
molecule. At first, the lower value was widely accepted as the most 
probable one (Herzberg, ref. 3 ) . Krieger and White (ref. 4) and 
Hirschf elder and Curtiss (ref. 5 ) have published tables of thermodynamic 
properties of high-temperature air based on this value. Gaydon (ref. 6 ) 
was perhaps one of the first advocates of the view that the higher value 
was the correct one. Subsequently a number of experiments were performed 
which confirmed Gaydon 1 s opinion, among them the measurements of strong 
shock waves in nitrogen made by Christian, Duff, and Yarger (ref. 7 ) and 
the detonation studies made by Kistiakowsky, Knight, and Malin (ref. 8 ). 
This rendered the work of references 4 and 5 obsolete, but shortly there- 
after Gilmore (ref. 9) computed the chemical composition, energy, entropy 
compressibility, and pressure of air as functions of temperature and den- 
sity based on the higher value for the dissociation of nitrogen. Later, 
Hilsenrath and Beckett (ref. 10) published a similar table of these prop- 
erties, but in much smaller increments of temperature and density. The 
calculations in both of these references (9 and 10 ) are highly refined 
in the sense that they not only account for the major components of air 
and their most significant energy states, but they also take into account 
a large number of the higher energy states which are infrequently excited 
even at high temperatures, and most of the very minor chemical components 
of air are included. Therefore these works are among the most detailed 
estimates for the thermodynamic properties of air which have been made. 
Also it is not likely that these works will become obsolete, since the 
values of all the important energy levels used in these calculations are 
now quite secure. However, it is desirable to have approximate expres- 
sions for these properties in closed form which can be solved without 
iteration. Such solutions would be particularly valuable, for example, 
in preparing tables to be used with the method of characteristics for 
calculating the flow of real air. 
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In contrast to the fairly satisfactory state of development in 
regard to thermodynamic properties, knowledge of the transport properties 
of air at high temperatures is in an elementary stage. It is generally 
agreed that an accurate calculation of the transport properties should 
be based on the rather rigorous theory of Chapman and Enskog for monatomic 
gases (ref. 11) . The extension of this theory given by Wang-Chang and 
Uhlenbeck (ref. 12) would be used for molecular gases with internal 
energy. For such calculations it is necessary to know the interaction 
potentials which exist between the gas particles so that the so-called 
collision integrals can be calculated. Hi rschf elder , Curtiss, and Bird 
(ref. 13) have developed the theory of intermolecular collisions to the 
stage where the collision integrals and the transport properties of inert 
molecular gases can be calculated with good accuracy. However, when the 
air dissociates, as it does at high temperatures, the atom-atom and atom- 
molecule potentials are needed, and these are not sufficiently well known 
to calculate the collision integrals. At the present time, calculation 
of the atom-atom and the atom-molecule potentials is being attempted by 
quantum mechanical methods. However, it may be some time before these 
solutions are available, and even when they are completed a formidable 
obstacle to the calculation of the transport properties remains. This 
occurs because two atoms may approach each other along any one of a num- 
ber of 'potentials, depending upon the principal quantum number of the 
electrons and the spatial orientation of the electron spin and orbital 
momentum vectors at the time of collision. These potentials are quite 
different, some are partly attractive and others are entirely repulsive. 

The solutions for all of these potentials are implicit in the quantum 
mechanical description of the problem. It is hoped that these potentials 
might be weighted by their probabilities to yield single effective poten- 
tials for each species, for if the atoms which interact along each dif- 
ferent potential must be treated as a separate species in the gas, the 
calculations of the transport properties will be laborious indeed, even 
with electronic computing. Problems of this type are further compounded 
when temperatures are considered where the air begins to ionize. In any 
event, there is an urgent need for an estimate of the properties of air, 
due to the demands created by the expanding realm of very high-speed 
flight. Therefore, an engineering approximation for the transport prop- 
erties of high-temperature air would be valuable in the interim while 
more exact solutions are being prepared. Even after more precise solu- 
tions are available, an approximation giving the transport properties in 
closed form will be just as desirable as it is for the thermodynamic 
properties of air. 

In view of the needs outlined above, it is the purpose of this paper 
to develop approximate expressions for the properties of air over the 
range of temperatures and pressures encountered by vehicles traversing 
the atmosphere at speeds up to escape velocity. The principle which shall 
be used in deriving these expressions is that they shall be made as simple 
in form as possible. For the thermodynamic properties this is accomplished 
by keeping o nly those terms which are necessary to yield a final result 
within a few percent of the more exact solutions. The properties of air 
which will be evaluated are the compressibility (i.e., the correction to 
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the ideal gas equation of state) , the energy, enthalpy, entropy, specific 
heat at constant pressure and at constant density, the speed of sound, 
the viscosity, thermal conductivity, and the Prandtl number. All of 
these properties will be evaluated for equilibrium conditions. Because 
of the finite reaction rates, these values will doubtless need to be modi- 
fied for processes which involve changes in state which are rapid compared 
to the rate of approach to chemical equilibrium. Nonequilibrium effects 
will probably be encountered in very high-altitude flight because the 
approach to equilibrium is slow at the low pressures experienced there. 
However, there is experimental evidence that equilibrium is essentially 
realized in flow of dissociating air under conditions encountered in 
flight at moderate altitudes and speeds (ref. 14) . In addition, 

Hirschf elder (ref. 15 ) has argued that heat transfer in pure conduction 
processes will correspond to equilibrium values if the reaction rate in 
one direction is rapid. Therefore, the thermodynamic and transport prop- 
erties of air which are based on equilibrium conditions should apply 
directly to some practical problems as well as being a convenient reference 
for the nonequilibrium values. 


a 


a i,bi 


Ai,Bi 


C 

Ci 

Cp 


Cv 


D 

D ij 


SYMBOLS 


speed of sound (zero frequency) 

stoichiometric coefficients for components and Bi 
components of a chemical reaction 
Sutherland’s constant (eq. (58)) 

specific heat per mol at constant density for component 

specific heat per mol at constant pressure 

partial specific heat per mol at constant pressure. 



i 


specific heat per mol at constant density 

dissociation energy per molecule, also diffusion coefficient 

binary diffusion coefficient for molecules of type i into 
molecules of type j 

base of natural logarithms, also electron charge 
electron 
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E 

Ei 

Eo 

Sn 

h 

H 

Hi 

I 

Io 

J 

k 

k 0 

k n ,k* 

k r 

K c 


K t 


Le* 


Zn 

m 

Mi 

M 

Mq 


energy per mol, also electric field strength 

energy per mol of component i 

energy per mol at zero absolute temperature 

degeneracy of the ith state 

degeneracy of the nth electronic state 

Planck* s constant 

enthalpy per mol 

enthalpy per mol of component i 

molecular moment of inertia, also ionization energy per 
molecule 

resonance potential for ionization 
rotational quantum number 

Boltzmann constant, also thermal conductivity 

reference coefficient of thermal conductivity (eq. (77)) 

coefficient of thermal conductivity due to molecular 
collisions 

coefficient of thermal conductivity due to chemical reaction 
chemical equilibrium constant for concentration units 
chemical equilibrium constant for pressure units 

DpC * 

partial Lewis number, - 

Mk* 

logarithm to the base e 

mass of a gas particle 

molecular weight per mol of component i 

mean molecular weight per mol of a gas mixture 

molecular weight per mol for undissociated molecules 
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n 


n (Ai ) , 

n(Bi), 

N 


No 

N + 

N 2 

NO 

0 

0 + 


0 2 

p 

Po 

p(Ai) , 
p(Bi), 

Pr 


Pr « 


Q 

Qt 

Qr 

Q v 

Qe 

Qc 


vibrational quantum number, also electronic quantum number, 
also concentration in moles per unit volume 

concentration of components Ai, Bj_, ... in moles per unit 
volume 

nitrogen atom, also atoms in general 

Avagadro number, molecules per mol 

nitrogen plus ion, also plus ions in general 

nitrogen molecule 

nitric oxide 

oxygen atom 

oxygen plus ion 

oxygen molecule 

pressure 

reference pressure, 1 atmosphere 

partial pressure of components A-^, B^, ... 

Prandtl number, 

’ Mk 

O-n’T) 

partial Prandtl number, -£■ — 

Mk* 

total partition function 
translational partition function 
rotational partition function 
vibrational partition function 
electronic partition function 

total partition function for a standard state of unit 

concentration, — Q 
' RT 
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Qp 



r 


r e 

R 

S 


Si 

So 

S ij or l 

S(i-j) J 


T 


uo 

U 

X 


x i 

x(A^) . . . 


Z 


a 

3 
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total partition function for a standard state of unit 
pressure, pQ 

total partition functions for components A-^, Bj_, ... 
distance between atoms 

distance between atoms where the potential is a minimum 
universal gas constant, energy per mol deg 
entropy per mol 

entropy per mol of component i at 1 atmosphere pressure 
collision cross section for undissociated air molecules 


collision cross section for particle i with particle j 
absolute temperature 

mean molecular velocity for molecule type i 
mean molecular velocity for undissociated air molecules 
potential energy between gas particles 
mol fraction 

mol fraction of component i 
mol fraction of component Aj_ ... 


pMo Mo 

compressibility, — — or 

pRT M 

molecular symmetry number (equal 2 for homonuclear diatomic 
molecules) , also polarizability 

Morse function constant (eq. (6l)) 


ratio of specific heats, 



Cv 


fraction of molecules which are dissociated or of atoms 
which are ionized 
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e i 


€ n 


5 


no 

\> 


Pi 


Po 

a 


P 

p 

s 


i,0 

■t,r,v,e 


energy 

of 

the 

ith 

state 

energy 

of 

the 

nth 

electronic state 


dimensionless distance parameter, — - 1 

r e 

coefficient of viscosity 

reference coefficient of viscosity (eq. (67)) 
mean free path for molecule type i 
reference mean free path (eq. (71)) 
vibrational frequency 
density of molecule type i 
reference density (eq. (70)) 
collision diameter 

Subscripts 

partial derivative at constant pressure 

partial derivative at constant density 

partial derivative at constant entropy 

indices referring to molecules type i and j 

indices referring to the contribution of translational, 
rotational, vibrational, and electronic energy modes, 
respectively 

THERMODYNAMI C PROPERTIES 


As a preliminary, a brief review of some of the results of 
statistical mechanics will be given. This will include the definitions 
of the partition functions and will summarize those relations between 
these functions and the thermodynamic properties of gases which will be 
used in the approximations to follow. 


NACA TN 4150 


9 


Partition Functions 


All of the thermodynamic properties of a gas may he calculated from 
its partition function. Consequently, the first step in determining the 
properties of air is to calculate the partition functions for the 
components in air. The partition function may he defined as 


where €j_ is the energy of the ith state of the gas particle and gp 
is the degeneracy, that is, the number of states of the particle which 
have this same energy level. The energy may he due to the translational, 
rotational, or vibrational motion of the particle, or to the motion of 
the electrons within the particle. The temperatures being considered in 
this paper are in all cases low enough that the excited nuclear energy 
states may be disregarded. The usual assumption is made that no coupling 
exists between the different modes of energy. Then the partition function 
may be expressed as the product 


The factors on the right side of equation (2a) are, respectively, 
the partition functions associated with the translational, rotational, 
vibrational, and electronic energy levels of the gas particle. Each 
factor is determined independently by an equation of the same form as 
equation (l). By the methods of statistical mechanics it is found that 
for diatomic molecules these factors are: 



( 1 ) 


Q — Qf^r^v^e 


(2a) 



(2b) 


00 


h 2 J( J+i) 




J=o 


(2d) 


(2e) 


n=i 
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The notations used above are common ones and the derivation of the 
equations may be found in any standard text on statistical mechanics, 
such as reference l 6 or 17. For monatomic particles, which have no modes 
of rotational or vibrational energy, the rotational and vibrational par- 
tition functions take the value unity. The translational and electronic 
partition functions for such particles take the same form as equations ( 2 b) 
and ( 2 e) , respectively. 

Consistent with the stated purpose of this paper, only those 
exponential terms are included in the electronic partition function 
(eq. (2e)) for which the energy levels, e n > are less than six times kT 
at the maximum temperature considered (15,000° K) . Actually, the levels 
are so widely split in this range that the closest to this cutoff is the 
sixth state of the atomic nitrogen ion, just a little more than four 
times kT at 15,000° K. 

Table I presents the atomic and molecular constants which were used 
in calculating the partition functions. The molecular constants for 
rotation, vibration, dissociation energy, and electronic energy levels 
were taken from Herzberg (ref. 3 ). The 9.76 electron-volt value for 
nitrogen dissociation is used, and it is assumed that the rotational and 
vibrational constants for all excited electronic states are the same as 
for the ground state. The atomic energy levels are taken from Moore 
(ref. 18). The constants have been rounded off to the nearest 0.5 per- 
cent, and the second and third electronic energy levels of atomic nitro- 
gen have been combined since they lie within 0.5 percent of each other. 

The same treatment applies to the second and third electronic levels of 
the atomic oxygen positive ion. 


The functions which are to be used directly in the calculations to 
follow are the logarithms of the partition functions. From the partition 
function constants (table I) and equations ( 2 a) through ( 2 e), these 
functions become 


7 n Q(Ng) 


1 to T - 0.42 

2 



to 0(02) 


1 to T + 0.11 

2 



l n 



11390 

+ 2 e 



In p 


(3a) 


(3b) 
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228 _ 326 

to Q(0) = ^ In T + 0.50 + Zn|5+3e T +e T + 


22800 _ 486 00 

5 e ^ + e ^ 


- In p 


(3c) 


In Q(W) =2 i n t + 0.30 + to f 4 + lOe 


27700 41500 

“T” + 6 e” T 


Zn p 


(3d) 


38600 58200 


to Q(0 + ) = ^ to T + 0.50 + to ( 4 + 10 T + 6 e 


-top 


(3e) 


70.6 188 .9 


Zn Q(W*) = — Zn T + 0.30 + Zn ( 1 + 3 e T + 5e T + 


22000 47000 67900 


5 e x + e 


T + 5 e T 


I - Zn p 


(3f) 


Zn Q(e~) = ~ Zn T - 14.24 - Zn p 


(3g) 


where T is the temperature in degrees Kelvin and p is given in 
atmospheres . 

It will he surmised that those components of air which are not 
represented by partition functions above are to be neglected in the 
approximations which are to follow. The absence of a partition function 
for nitric oxide may be found surprising, but it will be seen later that, 
to the order of accuracy being considered here, the formation of nitric 
oxide may be neglected in computing thermodynamic properties over a wide 
range of pressure and temperature conditions, including those conditions 
which will generally be encountered in high-speed flight through the 
atmosphere. This occurs because WO has about the same thermodynamic 
properties as an average for nitrogen and oxygen molecules, and the WO 
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formation does not greatly influence the equilibrium mol fractions of 
molecules and of atoms. The possible influence of nitric oxide on the 
transport properties of air will be considered later. 


Energy, Enthalpy, and Specific Heats for the 
Components of Air 


According to statistical mechanics, the energy and enthalpy per mol 
of pure gas are given by the following relations 


E-Eq _ [j] 1 


to Q\ 

_ T d to Qc 

(4) 

RT 

V 

8t / p 

dT 


H “ E o _ T 


to Q\ 

d Zn Qp 
= T p 

(5) 

RT 

V 

3T A 

dT 


The quantities Q c and Qp are the partition functions for the standard 
states of unit concentration and of unit pressure, respectively. These 
are related to the total partition function by 



(6a) 

Qp = PQ 

(6b) 


and they are functions only of temperature so that it is their total 
derivatives which are related to the energy and enthalpy as given in 
equations (k) and ( 5 ). The quantity E 0 is a constant representing the 
energy of the gas at zero absolute temperature. The choice of this level 
is arbitrary, but by convention E 0 is taken as zero for the molecules 
of nitrogen and oxygen. Then E 0 for the neutral atoms is just one half 
the energy of dissociation per mol of diatomic molecules , and E 0 for 
the ionized atoms is the sum of this dissociation energy and the energy 
of ionization. With all the ionization energy attributed to the ionized 
atoms, E 0 for the electrons must be taken as zero, of course. 

By definition the specific heats per mol of pure gas are 



(7) 


P 


(8) 
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From equations (2a), (2b), (2e), (4), and ( 5 ) it is seen that the 
energy and enthalpy per mol of gas due to translation and electronic 
excitation are given by 



and the specific heats given by equations ( 7 ) and. (8) can be expressed 



Equations ( 9 ) , (10), (ll), and (12) give that part of the energy, 
enthalpy, and the specific heats which is due to the translational and 
electronic energy for either atoms or molecules. The contributions of 
rotational and vibrational energy must also be included for the molecular 
case, of course. According to equations (2c) and (2d) the expression 



should be added to equations ( 9 ) and (10) in order to obtain the total 
energy and enthalpy for diatomic molecules, while 
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should, be added to equations (ll) and (12) to get the total specific 
heats for these particles. 

The entropy of a pure gas is related to its partition function by 

S 
R 

and from equations (5) and 
may be expressed 


= In Q + T 


Zn Q \ 

V ST I 


(15) 


(6b) , the entropy of the gas at unit pressure 


| = In Q p + 


H-E 0 

RT 


(16) 


The calculated values of the partition functions for a standard 
state of 1 atmosphere pressure, the energy, and the specific heat at 
constant density for the major components of air are given in table II 
in 500° K increments of temperature. The concentration standardized 
partition functions, the enthalpy, entropy, and the specific heat at 
constant pressure are not listed since they may be easily found from the 
properties tabulated and equations (6), (10), (12), and (l6). The tabu- 
lated properties for each component are given for the range in tempera- 
ture in which the component exists in air, at pressures between 0.0001 
and 100 atmospheres, up to the maximum temperature considered, 15,000° K. 
We shall now examine how these values may be used in computing the 
thermodynamic properties of air. 


Equilibrium Constants 


In order to determine the equilibrium mol fractions for the components 
of air, it will be necessary to calculate the equilibrium constants for 
the chemical reactions which occur. These chemical reactions may be 
expressed in the general form 


I 


aqAi 



(17) 


where the Aq are the reactants, the Bq the products, and aq and bq 
are their respective stoichiometric coefficients. The pressure equilib- 
rium constant for this reaction is defined in terms of the partial 
pressures 
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K t 


np bi (Bi) 

n P ai ( Ai ) 


( 18 ) 


and it is related to the partition functions by (see ref. 16) 

k p = " w + y b i zn aj - ^ q p^ a i^ 


where 


AEo =y biEo(Bi)- y a i Eo(A i ) 


(19) 


is the zero point energy of the products less the zero point energy of 
the reactants, both referred to their standard states. The reactions 
considered here are the dissociation of molecular oxygen and of molecular 
nitrogen and the ionization of atomic oxygen and of atomic nitrogen. 
According to equation (19) and the partition function constants given in 
table I, the equilibrium constants for these reactions are 


Zn Kp(C>2 — >2 0) 
Zn K p (N 2 -> 2N) 
Zn Kp(0 — ^ 0 + +e ) 
Zn K p (N li^+e") 


59,000 , 
T 

2 Zn Qp(0) - 

Zn 0p(0 2 ) 


(20a) 

113,200 

T 

■ 2 Zn Qp(N) 

- Zn Q p (N 2 ) 


(20b) 

158,000 ( 
T 

to Q p (0 + ) + 

to Qp(e") - 

Zn Qp(0) 

(20c) 

168,800 ( 

to Qp ( IP 1 " ) + 

to Qp(e-) - 

Zn Q p (N) 


T 

Jr 



(20d) 


The concentration equilibrium constant is defined by 


nn bl (Bj) 

nn ai (Ai) 


( 21 ) 


where n(Ap) and n(Bq) are, respectively, the concentrations of the 
chemical reactants and products. This quantity will also be needed for 
subsequent calculations, and it is obtained by replacing the pressure 
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standardized partition functions, Q p , with the corresponding concentration 
standardized partition functions, Q c (eq. ( 19 )). From equation (6) it 
is seen that 


K c 


= K p (RT) 


Ea^-Ebi 


(22) 


The logarithmic derivatives of the equilibrium constants will also 
be required later. From equations ( 5 ), (10), and- (19) these become 


d Zn K c _ AE 0 V" 1 /E-EqN f E -Eq\ 

~ " ~ + L\ RT ) “ L\ RT ) 

Bi A i 


( 23 ) 


T 


d Zn Kp 
dT 


d In K c V" 1 

+ ) in 

dT L 



(24) 


The equilibrium constants and their logarithmic derivatives for the 
reactions represented by equations (20a) through (20d) are listed as 
functions of temperature in table III. These quantities will now be 
used in calculating the component mol fractions and their derivatives. 


Calculation of the Equilibrium Mol Fractions 
and Their Derivatives 


The possibility that approximate solutions in closed form could be 
obtained for the properties of air suggests itself upon examination of 
the results of Gilmore (ref. 9)- His tables of the composition of air 
show that there are four chemical reactions of major importance. These 


are the dissociation of molecular oxygen and of molecular 
the ionization of atomic oxygen and of atomic nitrogen. 

nitrogen, and 

0 2 ~ ^ 2 0 

( 25 a) 

W 2 2N 

( 25 b) 

0 0+ + e" 

( 25 c) 

N If 1- + e" 

( 25 d) 


With one exception, all other reactions which occur yield component 
concentrations which are the order of 0.1 percent, or less, of the major 
components given by the reactions above. The exception is the formation 
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of nitric oxide, NO, which at sea level density may become as much as 
10 percent of the air around 5000° K (see ref. 9). However, even this 
much nitric oxide does not strongly influence the resulting thermodynamic 
properties of air, and at densities less than 0.01 normal sea level 
density, where the NO is less than 1 percent, the effects are very 
small. 

Two distinctive features of the chemical reactions given above are 
observable from Gilmore’s results. The first is that at all pressures 
the dissociation of oxygen is essentially complete before the dissocia- 
tion of nitrogen begins. This means that these two reactions can be 
treated independently for the purposes of our approximation. The second 
feature is that nitrogen and oxygen atoms ionize at about the same tem- 
perature and with about the same energy changes. (Note, e.g., the approx- 
imate equality of the equilibrium constants for the two ionization 
reactions, table III.) Consequently, it is possible to assume that once 
air is completely dissociated, all atoms constitute a single species 
which has the population weighted average properties of the nitrogen and 
oxygen atoms. Then the problem conveniently divides itself into three 
cases, and in each case only three components of the air need to be con- 
sidered simultaneously. These three cases are described as follows, 
where the ratio of nitrogen to oxygen has been taken as 4 to 1. 

Case I.- At relatively low temperatures, the three components of 
air which may exist simultaneously are molecular nitrogen, molecular 
oxygen, and atomic oxygen. If e is the fraction of the initial number 
of moles per unit volume which have become dissociated, then the partial 
pressures for these three components may be expressed 


p(N 2 ) = x(N 2 )p = M p 
1+e 

( 26 a) 

p(0 2 ) = x(0 2 )p = 2i£ l£ p 

1+e 

( 26 b) 

p(0) = x(0)p = — p 
1+e 

( 26 c) 


where the x(A^) are the mol fractions of component Aj_ . It is assumed 
that the gas components each behave like an ideal gas. Then the equation 
of state is 


p _ ZRT 
P ” Mo 


( 27 ) 


where the compressibility Z is 1 + e in this case. This quantity Z 
is the total number of moles per initial mol of undissociated air. It 
is also equal to the ratio of the initial molecular weight of 
undissociated air to the mean molecular weight, Mq/R. 
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The equilibrium constant for the oxygen dissociation reaction is 


K„(0 S 


20) is fp 

p(0 2 ) (l+e)(0.2-e) 


( 28 ) 


Then the fraction e is found by solving the quadratic equation (28) 



The condition of no dissociation which occurs at low temperatures is, 
of course, just the limit of case I where e is zero. 

Case II.- At intermediate temperatures the oxygen is completely 
dissociated and the nitrogen dissociation commences. For this case, e 
will be defined as the fraction of the initial air which dissociates into 
nitrogen atoms. Then the compressibility Z is 1.2 + e for this case, 
and the partial pressures of the components are 


p(N 2 ) 

= x(n 2 ) p = p 

1.2+e 

(30a) 

P(N) 

= x(N)p = P 

1.2+e 

(30b) 

p(o) 

= x(o)p = p 

1.2+e 

(30c) 

The equilibrium constant for 

the nitrogen dissociation reaction 

is 

K p (N 2 -» 2N) 

p 2 (N) 4e 2 p 

’ p(N 2 ) ~ (1.2+e) (0.8-e) 

(31) 

whence e is given by 



-0.4 + 

f — 


(32) 
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Case III.- At high temperatures, the dissociation of nitrogen is 
also complete, and the ionization of the atoms begins. For this case 
we assume that the atoms are a single homogeneous species symbolized 
by N. The term e will now be the fraction of the atoms which have 
been ionized, and then the partial pressures become 


P(N) 

= x(N)p = p 

(33a) 

1+e 


p(N+) 

= x(N+)p = P 

l+€ 

(334)) 

p( e_ ) 

= x(e")p = -i- p 
1+e 

(33c) 


Note that the compressibility Z in this case is 2(l+e). The ionization 
reaction equilibrium constant is 

Kp(M NUe-) = p(Mt) / P ' (e ' ) - (34) 

p p(N) 1-e 2 

and then e becomes 



(35) 


The derivatives of the mol fraction compressibility products, Zxj_, 
will be needed in addition to the mol fractions. For case I, from 
equations (28), (26a), (26b), and (26c) it is seen that 
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where the subscript p refers to the fact that the partial derivatives 
are taken at constant pressure. In order to find these same partial 
derivatives at constant density, the equilibrium constant is expressed 
in concentration units 


whence 


K c (0 2 20) 


4e g n 4e g p 

(l+e)(0.2-e) 0.2-e Mq 


d Zn K c m (g + 1 A /de\ 
dT \ e 0.2-e) W/p 


(40) 

(41) 


anfl the mol fraction derivatives at constant density take the same form 
as equations (37), (38), and (39) except that the quantity (8e/dT) p is 

replaced by the quantity (be /bn) p found in equation (4l). For case II 
a similar set of equations for the derivatives occurs 


d Zn k p = (2 _ 1 1 N fb£\ 

dT \e 1.2+e 0.8-e/ \bn) 

d Kc _ (2 1 N (be\ 

dT V € 0.8-e/ \8T/ p 

8Zx(M) o _ 2 8Zx(N 2 ) = 2 be 

8t 8t bn 

8Zx(o) = 0 

bn 


Finally, for case III, the equations are 


- Zn h + jJ)fbe) 

dT \ e 1+e 1-e/ w/ p 

d to K c /2 + j/\ /WN 

dT \e 1-e/ \8T/ p 

8Zx(N + ) _ 8Zx(e~) _ 8 Zx(N) ^ 8e 

8T ~ dT bH bl 


(42a) 

(42b) 

(43) 

(44) 

(45a) 

(45b) 

(46) 


where, as "before, the partial derivatives are given for constant pressure 
or for constant density depending on the corresponding partial derivative 
of e which is used. With these quantities in hand we are in position 
to calculate the energy, specific heat, entropy, and the speed of sound 
for air. 
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Energy, Specific Heat, Entropy, and Speed of Sound 
for Air in Equilibrium 


The energy per mol of air is simply the sum 


E = 



( 47 ) 


where the El are the energy per mol for component i. Generally in 
aerodynamic problems the energy per fixed mass of gas is needed rather 
than the energy per mol. The mass of gas in a mol of undissociated air 
(29 gm) will be used as the constant reference quantity for this purpose. 
Then, in dimensionless form, the energy per mol of initially undissociated 
air is 



( 48 ) 


while the dimensionless enthalpy per initial mol of air becomes 


— = — + Z 
RT “ RT + 

The compressibility, Z, and the dimensionless energy, ZE/RT, which 
have been calculated from the preceding equations, are listed in 
tables IV(a) and IV(b) and are graphed as functions of temperature in 
figures 1 and 2. The enthalpy is easily obtained from these values via 
equation ( 49 ). At this point it is desirable to compare the calculations 
of references 9 and. 10 with the results obtained above. It is not possi- 
ble to compare the numerical results directly since both references 9 
and 10 present the properties of air as functions of temperature at con- 
stant density rather than at constant pressure as is done here. However, 
the results of these two references are the same, and reference 10 tabu- 
lates the properties in close enough intervals so that constant pressure 
values can be obtained fairly accurately by interpolation. The values 
for compressibility and energy interpolated from reference 10 agree with 
the approximate solutions within 5 percent in a.ll cases, and generally 
within 2 percent, so that on the scale of figures 1 and 2 the differences 
are hardly discernible. In view of this agreement, it is concluded that 
the solutions presented here are useful approximations for the thermody- 
namic properties of air. Therefore we shall proceed to use those approx- 
imations to calculate the entropy, specific heats, and speed of sound in 
air. 



The entropy per initial mol of air is simply obtained from the 
entropies of the components of air through the relation 
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where p Q is the reference pressure for the standard state, in this case, 
1 atmosphere. The entropy values are listed in table IV(c) and are shown 
graphically in figure 3 • 

It is necessary to take the derivatives of the products ZE and ZH 
for specific heats if the latter are to have their usual meanings, that 
is, the change in energy which occurs in a fixed mass of gas per degree 
temperature change. Thus the constant density specific heat per initial 
mol of air is given by 


ZC V 

R 


= I fdZE\ 
R VST A 



(51) 


where Cp is the derivative of energy for component i, that is, dEi/dT. 
The corresponding specific heat for constant pressure is given by 


ZC P _ 1 f] SzhV 
R RVST/p 



(52) 


The specific heats calculated from equations (51) and (52) are 
listed in tables IV(d) and IV(e). It may be noted that when chemical 
reactions occur, the difference between the specific heats per mol is 
not equal to the gas constant, as in an ideal gas, but is given by 



Both of the specific heat functions have the same general features. 
These features are illustrated in figure 4 where the specific heat at 
constant density is given as a function of temperature for pressures 
of 0.0001, 0.01, and 1 atmosphere. At low temperature C v increases 
from 5/2 R to 7/2 R as the vibrational modes of energy become excited. 
Then, with increasing temperature, the specific heat goes through three 
distinct maxima where the chemical components change most rapidly with 
temperature; the first maxim-urn is due to the oxygen dissociation reaction, 
the second to the nitrogen dissociation, and the third to the ionization 
reactions. As pressure decreases, these maxima increase in sharpness and 
in magnitude and they shift to lower temperatures. It has been pointed 
out (ref. 16) that the values determined from the partition functions 
cannot be highly accurate because the second derivatives of the partition 
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functions are required, and the error in neglecting higher order terms 
is magnified in the process. However, the curves of figure 4 would 
merely be shifted slightly in temperature and have somewhat different 
peak values if the calculations were more exact. Certainly the qualita- 
tive features of the specific heats should be correctly given. 

The specific heat values determined above enable us to calculate 
the speed of sound in air. This speed of sound will be defined by 


a 2 = 



( 5 ^) 


that is, the limiting value as the frequency of sound approaches zero. 
The partial derivative given in equation (54) will not be calculated 
directly, since it is not convenient to treat the entropy as an inde- 
pendent variable which can be held constant (see eq. (50))- However, 
equation (54) may be transformed with perfect generality into the form 


a2 = 7 W T 


Bp 


(55) 


where 7 is the ratio of the specific heats, Cp/C v . Equation (55) 
in turn, equivalent to 


a 2 = -7 


Bp 

St 


(56) 



which, from the equation of state (eq. (27)), may be expressed in terms 
of variables which have already been calculated 



(57) 


The dimensionless speed of sound parameter, a 2 p/p, is listed in table IV(f) 
and is plotted as a function of temperature in figure 5* The second term 
on the right side of equation (57) is generally near unity, so that fig- 
ure 5 is also indicative of the variation in 7 with temperature. 

The range of pressure and temperature in which the dissociation of 
oxygen, the dissociation of nitrogen, and the ionization of atoms are the 
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important chemical reactions is indicated by dividing table IV into the 
regions case I, case II, and case III, respectively. Figure 6 presents 
the same information in a graphical form from which it is convenient to 
judge the chemical state of air for a given pressure and temperature. 

It will be noted that boundaries between the chemical reactions corre- 
spond to end points in the characteristic cycles that occur in the 
thermodynamic parameters graphed in figures 1 through 5* 


Aerodynamic Considerations 


The thermodynamic properties obtained at this point are those 
required to perform calculations of inviscid air-flow problems by the 
method of characteristics. These properties are given for a range of 
temperature from 500° to 15,000° K and of pressure from 0.0001 to 
100 atmospheres. It is of interest now to examine the altitude and 
velocity at which these conditions will occur in flight. A grid of the 
pressure and temperature at the stagnation point of a body in flight is 
shown in figure 7 as a function of flight altitude and velocity. The 
stagnation enthalpy per unit mass was simply taken as one half the veloc- 
ity squared, and the stagnation pressure was related to the static pres- 
sure (and thus to altitude) with the results of Feldman (ref. 19) , who 
computed the pressure ratio developed across normal shock waves in air 
at various altitudes. Generally lower temperatures and pressures will 
be attained at regions other than the stagnation region, so the range of 
these variables will be adequate for such cases also. At very low pres- 
sures, such as occur at high altitudes or in expanded flow regions, aero- 
dynamic effects may often be disregarded. However, if necessary, the 
approximations of this report may be extended to pressures below 0.0001 
atmosphere. It is not advisable to extend them to pressures much greater 
than 100 atmospheres, since there occurs an increasing overlap of the 
chemical reactions in air as pressure increases and, in addition, the 
nitric oxide becomes such a large component of the air that it must be 
considered. As a result of this overlap in the chemical reactions, the 
calciilated properties of air show irregularities in the region of transi- 
tion from one reaction to another. Slight irregularities can be observed 
at the transition regions for 10 and 100 atmosphere pressures (figs. 1 
through 4), for example. However, the deviations are small at these pres- 
sures and it can be seen from figure 7 that the thermodynamic properties 
of air can be closely approximated in closed form over the range of con- 
ditions of current interest in aerodynamics. Some of the results of these 
approximations will next be used in estimating the transport properties 
of high-temperature air. 
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TRANSPORT PROPERTIES 
Collision Cross Sections 


Consider first some qualitative aspects of the collisions between gas 
particles. The particle trajectories are influenced by a potential U 
which is negligible at long range , which may be attractive or repulsive 
at intermediate range, but which always becomes strongly repulsive at 
very short range. A particle with kinetic energy kT will not be greatly 
deflected if it passes only through that part of the potential field 
where |u | « kT. On the other hand, the particle will be deflected appre- 
ciably if it passes where |u | » kT. The direction of the deflection is 
unimportant so far as transport properties are concerned, for it is the 
absolute value of the deflection which determines the change in mass, 
momentum, and energy fluxes caused by the collision. To a first approxi- 
mation, the absolute value of this deflection is independent of the sign 
of the potential and the effective collision diameter a is the order 
of the largest distance where U = ±kT. The effective collision cross 
section S will be defined as jrcr 2 . 

In the rigorous treatment of the transport properties of gases, the 
effective collision cross section S is found to be an integral function 
of the deflection angles produced by a collision, and this integral is a 
function of the relative velocity of the colliding particles. The 
so-called "collision integral” is a function of temperature only which 
is S times a velocity function integrated over all velocities. Thus 
the collision integral may be thought of as an average effective colli- 
sion cross section, and the transport coefficients can be related directly 
to these integrals. However, not all the interparticle potentials have 
been developed which are needed to calculate the collision integrals for 
air. In the present paper then, plausible estimates of the effective 
collision cross sections will be used to determine the mean free paths 
for hard elastic spheres, and for such particles the transport coeffi- 
cients can be related to these mean free paths (ref. 1). The effects of 
the interaction potentials will be taken into account by letting the 
spherical cross sections be appropriate functions of temperature. 

The effective cross sections for collisions between diatomic 
molecules can be obtained quite accurately by the collision integral 
method. However, at high temperatures, the very steep repulsive portion 
of the intermolecular potential is penetrated so that the molecules do 
behave essentially like hard spheres. Consistent with the approximations 
which are to follow, it will be sufficient to use the Sutherland formula 
for the molecular cross section S 0 


where C is 112° K and S^, the molecular cross section at infinite 
temperature, is 3-l4xl0" 15 cm 2 . 


(58) 
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two 

U(r) 


0 


U Q (a) = -kT (59) 

and these diameters will he used to evaluate the coefficients of momentum 
Anri energy transfer. A somewhat deeper penetration of the potential is 
normally required for a collision to affect the particle flux, so that 
the diameters a* which will be used to evaluate the diffusion 
coefficients are assumed to be given by 

Uo(a’) = -2kT (60) 

It may be pointed out that Hirschfelder and Eliason (ref. 20) have 
examined the relation between values of the transport coefficients given 
by the hard sphere model and by the more rigorous collision integral 
method. They find that U(cr) and U(ff') are about -0.6kT for a wide vari- 
ety of attractive potentials, that U(cr) is about 0.9kT and U(cr’) is 
about 1.6kT for a similar variety of repulsive potentials. If all the 
interparticle potentials were known, it would be simple to use these 
criteria to obtain a weighted average collision diameter. All the poten- 
tials are known for two normal hydrogen atoms (ref. 21) and it is found 
for this case that equations (59) and (60) yield values for cr and a’ 
which agree with the weighted average collision diameters within 8 per- 
cent over the range of temperatures from 1000° to 15,000° K. Of course, 
there is no assurance that these same relations will hold as closely for 
collisions between the atoms in air. In fact the average collision diam- 
eters for normal oxygen and nitrogen atoms will probably be overestimated 
by equations (59) and (60) , since the shallow intermediate potentials 
(such as Ui and U 2 , sketch (a)) must be considered for these atoms, 
whereas they do not occur for hydrogen. The effect of these intermediate 
potentials will be partly compensated for by the fact that some of the 
atoms will be in excited electronic states which have collision diameters 
the order of three times larger than the normal atoms (ref. 20). The 


For atomic collisions, the picture is complicated by the fact that 
atoms may approach each other along any one of a number of potentials. 

For example, the potentials between two nor- 
mal nitrogen atoms are shown qualitatively 
in sketch (a) (see ref. 3) • The lowest lying 
of these potentials, designated Uq, has the 
lowest total electron spin and it is the one 
normally responsible for the vibrational 
energy levels observed in the stable molec- 
ular state. Therefore U Q can be expressed 
quantitatively from experimental spectro- 
scopic data. Unfortunately, the higher lying 
potentials for the atoms in air are not known 
quantitatively at present, so we are reduced 
to estimating an average collision cross 
section for all of the potentials by means 
of the known lowest lying potential Uo. 

For this purpose it is assumed that the 
Sketch (a) collision diameters cr are given by 
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fraction of atoms in excited states is small over most of the temperature 
and pressure range considered so that collisions between excited parti- 
cles are relatively rare . The encounters between an excited and a normal 
atom are the ones which significantly influence the mean free paths, and 
the cross sections for such collisions are about four times normal. The 
fraction of each constituent in air which is excited is listed as a 
function of temperature in the following table: 


T, 

°K 

Fraction of 
oxygen atoms 
excited 

Fraction of 
nitrogen atoms 
excited 

Fraction of 
ionized atoms 
excited 

4,000 

0.002 

0.002 


6,000 

.013 

.025 

0.012 

8,000 

.032 

.081 

.032 

10,000 

.055 

.153 

.060 

12,000 

.079 

.228 

.089 

14,000 



.120 


The first two excited states of atomic oxygen and of the atomic nitrogen 
plus ion are very close to the ground state (see table I) and so they 
have not been counted as excited states. It can be seen that the frac- 
tion of excited particles in air will generally be less than 10 percent 
and in the very worst case considered about 20 percent of the atoms will 
be excited. This occurs at 100 atmospheres pressure where the ionization 
reaction is repressed until temperatures beyond 12,000° K are reached and 
23 percent of the nitrogen and 8 percent of the oxygen atoms are excited. 
For a pure gas in which excited particles have three times normal colli- 
sion diameters, the average mean free path is decreased 28 percent when 
10 percent of the particles are excited, and is decreased 45 percent 
when 20 percent of the particles are excited. Thus the error introduced 
by neglecting the long-range forces of the excited atoms is probably the 
same order of magnitude and opposite in effect to the uncertainties caused 
by neglecting the shortened range of the shallow, intermediate collision 
potentials . 

The lowest lying atom-atom potentials may be approximated by the 
Morse function 



where D is the dissociation energy and (; is the dimensionless distance 
parameter (r/r e -l)j r e is the interatomic distance at which U Q is a 
min imum (sketch (a)); (3 is a constant related to the vibrational fre- 
quency v of the stable molecule and is given by v(jt 2 l/2D) 1/2 . Although 
the Morse function is not very accurate at long ranges, it will be needed 
only at rather high temperatures where the potential is fairly well 
described for |u D I > kT. The Morse function constants were taken from 
Herzberg (ref. 3) snd these are listed in the following table: 
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Atomic 

particles 

D/lc, 

6 K 

3 

angstroms 

0-0 

59,000 

3.24 

1.207 

N-N 

113,200 

2.96 

1.094 

N -0 

75 , 4 oo 

3- 18 

1.151 

0-0 + 

75,200 

3- 18 

1.123 

N-N+ 

101,200 

2.94 

1.126 


From equations ( 59 ) and. (6l) , the momentum and energy transfer 
collision diameter a is given by 

^ = 1 - 3 to (l-/TTf) (62a) 

and by equations (60) and (6l) the diffusion collision diameter a’ is 
given by 


£l = 1 - S to ( 1 . - ji-2f) 

The atom-molecule collision diameters will he taken as the arithmetic 
average of the atom-atom and of the molecule-molecule collision diameters. 
This assumption corresponds to the concept that the collision diameter is 
a measure of the effective range of the electron distribution about the 
molecule or atom, and that a collision occurs whenever these electron 
distributions overlap. 

With the onset of ionization, several additional types of collisions 
must be considered. These are the ion-atom, electron-atom, and the ion- 
ion, ion-electron, or electron-electron collisions. The ion-atom colli- 
sions may be treated in the same manner as the atom-atom collisions, since 
the spectroscopic data for ionized molecules are available to be used in 
setting up the Morse potential function between the ion and the neutral 
atom. The Morse function constants for these cases, which occur in air, 
are also listed in the preceding table. The ion-ion, ion-electron, and 
electron-electron collision diameters all depend on long range electro- 
static forces of identical magnitude. The criterion for the collision 
diameter will be taken as before |u ( cr) | = kT. The function U in this 
case is just the coulomb potential between two charges, and the collision 
diameter thus becomes 



— 

kT 


a 


(63) 
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For electron- atom collisions, the atom will he polarized by the 
approach of the electron and there will result a charge -dipole type inter- 
action. In order to calculate the magnitude of this interaction, the 
polarizability of the atom, a, will he needed. The polarizability is 
defined as the dipole moment produced in the atom hy an electric field, E, 
of unit strength. The interaction energy of the induced dipole moment 
with the field # is aE 2 /2, and when the field due to the presence of the 
electron is substituted for E, this energy becomes 


U = 


e 2 a 

2r 4 


( 64 ) 


It will be noted that the polarizability need not be known with 
great precision, since the collision diameter depends only on the fourth 
root of this quantity. Experimental values of polarizability are only 
given for the molecular state of oxygen and nitrogen. However, Joos 
(ref. 22) gives an approximate method of calculating the polarizability 
of atoms in alternating fields which, extrapolated to steady state fields, 
becomes 


a 



(65) 


where m in this case is the mass of an electron and I Q is the 
resonance potential for ionization of the atom. This resonance potential 
is 9-11 e.v. for oxygen atoms and 10.28 e.v. for nitrogen atoms (ref. 23). 
Using these values, one finds the polarizability is 13.2xl0~ 25 cm? and 
10. 3xl0" 25 cm? for oxygen and nitrogen atoms, respectively. 

The collision diameters will enter the calculations to follow in 
the form of a ratio with the collision diameter for two diatomic molecules 
at the same temperature. These ratios were calculated by the methods 
outlined above for those collisions which occur between the major com- 
ponents in high-temperature air. The results are given in table V as 
functions of temperature. Wherever two or more types of collision give 
the same result within 5 percent, there is not much point to considering 
them separately, in view of the other approximations involved. There- 
fore, for such cases, the values have been weighted and averaged, and 
only those final averages are given in table V. The atom-molecule cross 
sections S(0 2 -0), S(N 2 -N) , and S(N 2 -0) are very close and these are all 
designated by the single symbol S(N 2 -N) . Similarly the atom-atom and 
atom-ion cross sections S(0-0) , S(N-N), S(N-O) , S(N-E + ), and S(0-0 + ) 
are all grouped under the notation S(N-N) . The atom-electron cross 
sections S(N-e) and S(O-e) are given the single average value listed 
under S(N-e) . In the same way, all the atom-molecule and atom-atom 
diffusion cross sections are given by S T (N 2 -N) and S f (N-W), respectively. 
In view of the order of the approximations involved, the electron-atom 
cross sections for diffusion are not differentiated from the cross 
sections for momentum and energy transport. 
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Viscosity 


The viscosity of air will he calculated from the simple summation 
formula for a mixture of hard spherical molecules 



Pi u i^i 


( 66 ) 


In view of the other approximations involved, a more sophisticated 
formulation than this does not seem warranted. The constant 5«/3 2 is 
chosen to agree with the value which accounts for the persistence in 
velocity of the higher speed molecules (ref. 1). It will be convenient 
to use, as a reference value, the viscosity which air would have at the 
same temperature if the molecules did not dissociate or ionize 


4° - ^ Po^o 


(67) 


In cgs units, this reference viscosity becomes 


% 


1.462x10“ 5 


nTt 


gm 


1+112/T cm- sec 


( 68 ) 


where T is the absolute temperature in degrees Kelvin. The ratio of 
mean velocities is 


u i 


uo 



(69) 


and the ratio of densities, where the total number of moles per unit 
volume is the same, is 


% 

Mo 


Xi 


(70) 


Under these same conditions, that is, the same mol concentration, the 
ratio of mean free paths is (ref. 1) 
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Then the ratio of viscosity to the reference viscosity becomes 



i 




(72) 


Before equation (71) can be applied to the case where ionization 
occurs, it must, be observed that the atom-electron collisions are not 
to be counted in evaluating the atom mean free paths, nor are the ion- 
electron collisions to be counted in evaluating the ion mean free paths. 
This is because the mean free path which occurs in the derivation of 
viscosity is the mean distance between those collisions which cause a 
relatively large change in momentum of the particle being considered, 
either in direction or in magnitude. Normally all collisions qualify 
in this sense, but when a heavy particle like an atom or an ion collides 
with an extremely light particle like an electron, there is practically 
no change in the momentum of the heavy particle, and this momentum is 
carried intact to the first collision with another heavy particle. 

A feature of the ionized gas which should be mentioned is that charged 
particles under the influence of coulomb potentials may be deflected by 
the cumulative effect of many long-range collisions. Then, as pointed out 
by Cohen, Spitzer, and Routly (ref. 24), the effective collision cross 
section is larger than for a single collision at short range. However, 
a charge -shielding effect limits the range of the coulomb potential and 
where charge density is large the correction to the collision cross sec- 
tion is small. For the range of variables treated here, collisions with 
neutral particles generally predominate when the charge density is small. 
Thus, in the first approximation, the mean free path is considered to be 
a consequence of single, independent collisions as in equation (7l)» 

The coefficients of viscosity for high-temperature air have been 
calculated from equations (68), (71), and (72) and the results are pre- 
sented in table VI (a) and in figure 8. The terms in the summation of 
equation (72) will also be used in the evaluation of the coefficient of 
thermal conductivity which follows. 


A theory for the thermal conductivity of a chemically reacting gas 
was perhaps first outlined by Nemst (ref. 25) • The form in which this 
theory has been developed by Hirschf elder (ref. 15 ) will be used here to 
calculate the conductivity of high-temperature air. In this theory the 
energy transfer through the gas is treated in two independent parts. The 
first part is the energy transferred by molecular collisions, and this 
mode of energy transfer is the one responsible for the ordinary thermal 
conductivity of nonreacting gases. The second part is the energy trans- 
ferred by diffusion of the molecular species and the reactions which 
occur as the gas tends to maintain itself in chemical equilibrium at each 


Thermal Conductivity 
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point. We shall first turn attention to the calculation of energy 
transferred ty molecular collisions. If one invokes Eucken' s assumption 
(ref. 2 6 ), that the internal energy is distributed among the gas particles 
independently of their velocity distribution, then the simple linearized 
expression for the coefficient of thermal conductivity in a mixed gas 
becomes 

• ^ = flZ PiUi W (2 Ct + Cint )i (73) 

i 

where Ct is the specific heat per mol due to translational energy 
and Ct n t the specific heat due to the internal energy. It is noted 
that 


Ci - (Ct + ^inth 


(74) 


and 



whence equation (73) becomes 



i 


(75) 

(76) 


It will be convenient to define a reference coefficient of thermal 
conductivity just as was done for the viscosity, that is, the value air 
would have if it did not vibrate, dissociate, or ionize. This coefficient 
will be 


. 19 R 

ko = T * "o 


or in cgs units 


ko = 1.364 r) 0 


joule 
cm- sec °K 


where T) 0 is in 


gm/cm-sec. 


(77) 


(78) 


The ratio of the thermal conductivity coefficient to the reference 
coefficient becomes 
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i + 19 -l 


( 79 ) 


The first factors in each term are just the terms already calculated for 
the viscosity. It may he noted that the same collisions are to he neg- 
lected in this case as for viscosity. This is because these electron- 
atom and electron-ion collisions have only a small effect on the magnitude 
and direction for the energy flux vector of the heavy particles. Actually, 
all of the electron collisions could have heen neglected for the purpose 
of computing viscosity, for the electrons carry a negligible fraction of 
the total momentum transferred, due to their small mass. However, because 
of its high velocity, an electron transports a large share of the kinetic 
energy, in fact much more than a heavier gas particle. Therefore the 
heat conduction terms for the electrons are needed in equation (79). 

The second mode of energy transfer, which takes place whenever the 
gas undergoes a chemical reaction, is due to the diffusion of the chemical 
species. These particles then react with one another, giving off or 
absorbing the heat of reaction and causing heat transfer which may be 
considerably larger than the ordinary heat transfer due to molecular 
collisions. Hirschfelder (ref. 15) has formulated this problem, but in 
terms of the multicomponent diffusion coefficients which are difficult 
to estimate. However, Butler and Brokaw (ref. 27) have shown how 
Hirschfelder' s results can be modified to make use of the binary diffu- 
sion coefficients instead, and their solutions are in convenient form 
for computational purposes. Moreover, Butler and Brokaw show that 
Hirschfelder' s method of predicting thermal conductivity agrees very well 
with the experimental results for a number of gases which dissociate or 
react chemically at ordinary temperatures. Therefore we shall be able 
to apply this method with some confidence. That part of the coefficient 
of thermal conductivity which is due to the chemical reaction will be 
designated k r , and according to the results of reference 27, this may 
be expressed as 


where Dpj is the binary diffusion coefficient for molecules of type i 
into molecules of type j, and the a± are the stoichiometric coefficients 
of components Aq in the chemical reaction written in the form 



( 80 ) 



i J 
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^a-jAi =0 (81) 

i 

According to reference 1, the binary diffusion coefficients are given by 


1 

HDij 



(82) 


Now the reference coefficient of thermal conductivity ( eq. ( 77) ) can be 
put in the form 


k = 22 R 
64 N 0 S 0 

Substituting equation (82) in equation (80) and dividing by equation (83), 
one obtains 



kp 

^o 



i <5 


12 f 

T 

95 V 


d In Kp^ 2 
dT ) 


MjMj 


LMo(M i+ Mj)J 


1/2 q’- • o 

s ij ai 


^ ~ ( a i x ,1~ a ,1 x i) 


J M 


(84) 


In order to simplify the computations somewhat, the differences in 
mass between oxygen and nitrogen atoms will be neglected. Then, since 
some of the collision cross sections may also be given the same value, 
the double sum m ations in the denominator of equation (84) take on the 
relatively simple form as follows : 

For case I (oxygen dissociation reaction only) 

S» (N 2 -N) f [x(0)+2x(0 2 ) ] 2 + 4x(N 2 )| + x(N 2 ) (85 - 

*'/3 S Q 1 x(0)x(0 2 ) x(0) J n/2x(0 2 ) 

i j 

for case II (nitrogen dissociation reaction only) 
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S« (N 2 -N) 

-s/3 So 


{ 


[x(N)+ 2 x(N g ) ] 2 + x( 0 ) 

x(W)x(W 2 ) x(N 2 ) 


} 


S* (N-N) 2x(0) 
S 0 x(N) 


( 86 ) 


and finally for case III (ionization reaction only) 



"l S» (N-N) 1 S' (W-e-) "[ [x(N)+x(H +)] 2 

_2 S 0 + 230 S J x(N)x(N + ) 


‘l S» (N-N) 
_2 S 0 


(87) 


i j 


The total coefficient of thermal conductivity is just the sum, k n +k r . 
The values of this sum were calculated for air and the results are pre- 
sented in table VI (b) and graphed in figure 9 . The coefficient of thermal 
conductivity has about the same functional features with respect to tem- 
perature and pressure that were observed for the specific heat (fig. 4) . 
The Prandtl numbers were calculated from the relation 


and they are listed in table VI (c) and are graphed as functions of 
temperature in figure 10 . 

As noted, nitrogen and oxygen have approximately equal collision 
cross sections, so that up to the point where ionization begins, air can 
be treated essentially as a two-component mixture of atoms and molecules. 
Then it is possible to characterize the system with a single diffusion 
coefficient, and the differential equations of fluid flow may be analyzed 
while keeping the terms describing heat transfer by chemical reaction 
separate from the terms describing ordinary thermal conduction. This 
method is followed, for example, by Fay and Riddell (ref. 28) in their 
analysis for the heat transfer to the stagnation region of blunt bodies 
in high-speed flight. In using such an analysis, it must be realized 
that the thermodynamic and transport coefficients are defined somewhat 
differently than in the present paper. The coefficients which appear 
in reference 28 will be called the "partial coefficients" and they will 
be designated here by a superscript prime. The partial specific heat, 
the partial coefficient of thermal conductivity, the partial Prandtl num- 
ber, and the partial Lewis number, respectively, are related to the 
quantities defined in this paper as follows: 



( 88 ) 
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J 

The partial Lewis number, Le 1 , appears as a coefficient for the 
chemical reaction terms in reference 28. The diffusion coefficient D 
is the ordinary binary diffusion coefficient given by equation (82) 

and from equation (83) a dimensionless group is 


. 0.309 

koMo |_ Sq 


(90) 


The Lewis number is conveniently calculated from the expression 


(91) 


Table VII lists the calculated values for the partial coefficients 
ZC p « /R, k n /k Q , Pr f , and Le f . The numbers Pr* and Le 1 are also graphed 
as functions of temperature in figure 11. The partial coefficients are 
not given for temperatures where ionization occurs because the air is 
then at least a three -component system (neutral atoms, ionized atoms, 
and electrons) and the partial coefficients, as defined, do not apply to 
this case. It can be seen that Pr 1 is a relatively constant quantity 
with an average value about 0.72. The partial Lewis number, on the other 
hand, decreases as dissociation proceeds to completion. From the result 
that the partial Prandtl number is approximately constant, it follows 
that the factor Cp T /Mk 1 decreases as rp 1 . The factor Dp/M increases 
at about the same rate (see eq. (82) or ( 90 )), and the product of these 
two factors is relatively constant; thus the partial Lewis number varies 
approximately as M, that is, it decreases by a factor of about 2, from 
about 1.4 to 0.6. 

A list of conversion factors is given in table VIII for convenience 
in converting the tabulated and graphed parameters into dimensioned units 


ZDpR ZCp'/R 
k 0 Mo Z,k n /kQ 
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DISCUSSION OF RESULTS 


It is convenient to summarize the physical processes which influence 
the transport properties of air in a discussion of the variations in 
Prandtl number with temperature. At low temperatures the air is like a 
pure diatomic gas with a specific heat about 7/2 R. From equation (88), 
the Prandtl number for such a gas is 14/19 or 0.74. As temperature 
increases, vibrational energy is excited in the molecules so that Cp 
approaches 9/2 R, while from equation (76) the thermal conductivity 
approaches (23/4) (Rt] 0 /Mo) . For this limit the Prandtl number becomes 
18/23 or O.78. 


At still higher 
temperatures the oxygen 
dissociates and both 
Cp and k go through 
the pronounced maxima 
shown in figures 4 
and 9* To a first 
approximation, k is 
about proportional 
to Cp just as for a 
nonreacting gas . How- 
ever, due to the influ- 
ence of the chemical 
reaction, the maximum 
in k occurs at 
slightly lower temper- 
atures than the maximum 
in specific heat as 
indicated in sketch (b) . 


1.0 — 



k/k f 


Temperature, T 

Sketch (b) 


The viscosity coefficient is not greatly influ- 
enced by the oxygen dissociation, so from equation (88) it is apparent 
that the Prandtl number will follow an S shaped function as shown in 
sketch (b) . This shape for the function can be seen in the curves of 
figure 10. 


At the temperature where oxygen dissociation is essentially complete 
and nitrogen dissociation has not yet begun, the Prandtl number can be 
calculated for a two -component mixture from the known conditions that 
the mol fractions are x x = 2/3 and x 2 = 1/3 and that the molecular 
weights are related by = 2M 2 (where subscript 1 refers to nitrogen 
molecules and subscript 2 to oxygen atoms) . If the collision cross 
sections are all assumed equal, the mean free paths are related by 


A2 

Ai 


4+ n /6 

2+\T3 


( 92 ) 
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Then the Prandtl number becomes 




23 R 

5 HI. 


1 + ^ 1 ^ 
4 Ax 

_R_ (23 15J2 Aa\ 

Mi \4 8 AiyL 


0. 72 


(93) 


which checks closely with the values for Prandtl number shown in figure 10 
at the transitions from oxygen to nitrogen dissociation. 

As nitrogen dissociation proceeds, the Prandtl number again follows 
an S shaped function of temperature for the same reasons which were 
outlined for oxygen dissociation (sketch (b)). The fully dissociated 
air is like a pure monatomic gas with Cp about 5/2 R and k equal 
15/4 (Rq/M) , so that the Prandtl number approaches 2/3. The Prandtl 
numbers of figure 10 go through this value of 2/3 at a temperature where 
the dissociation is essentially complete but the ionization is still 
negligible (Z = 2, fig. 1). 

Up to this point, the Prandtl numbers lie within the range from 0.6 
to 1.0 in agreement with the conclusions of reference 29. However, when 
ionization begins, the Prandtl number may drop to somewhat lower values. 
For very small degrees of ionization it can be shown that the ratio of 
electron to atom mean free path is about «/2 and that the coefficient of 
viscosity is influenced very little by the presence of the electrons. 
However, the thermal conductivity is greatly increased by the electrons 
because of their high thermal velocity, and the effect is illustrated in 
figure 9 by the abrupt change in slope for k where the ionization reac- 
tion first begins. If the fraction of ionization, e, is so small that 
chemical reaction terms may be neglected 


k = 15 JL 
4 Mj. 


n 



(94) 


where Mi is the atom mass and M 2 the electron mass. It can be seen 
that k increases rapidly with the onset of ionization because the 
factor (2M 1 /M 2 ) 1 ^ 2 is large, about 230. Accordingly the Prandtl number 
is found to vary as 



3(l+230e) 


(95) 
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A very small value of e is sufficient to reduce the Prandtl number 
considerably below 2/3, the value for a pure monatomic gas. As the 
ionization proceeds, however, the chemical reaction terms which were 
omitted in equations (94) and (95) become predominant such that Cp 
becomes approximately proportional to k again, and the Prandtl number 
levels out at about 0.3 as shown in figure 10. At this point another 
factor gradually predominates, namely the decrease in coefficient of 
viscosity which occurs as the mean free paths become very short as a 
result of the long range coulomb forces acting on the charged gas parti- 
cles. Consider, for example, the completely ionized gas. The mol frac- 
tion of ions (subscript 1) equals the mol fraction of electrons 
(subscript 2) and the collision cross sections are about equal. Then 
the ratio of mean free paths is 


^2 * s /"2 

M ~ I+J2 


(96) 


where, as before, the electron collisions do not count toward the ion 
mean free paths. The ratio k/t} is approximately given by 

(97) 

and with Cp/M about the Prandtl number becomes very small 


k „ 15 _R_ {2 /*£ 

4 ^ 1+ </2 J M 2 


£p J2 _ h 1+ -J 2 

M k " 3 n/2 


= 0.014 
M x 


(98) 


as shown in figure 10. With the assumption that the Prandtl number 
retains its usual significance, boundary-layer regions in highly ionized 
flow should be much better heat conductors than in the flow of neutral 
gas particles. This result is of theoretical interest even though fully 
ionized air flow will not be obtained at flight velocities below escape 
speed (figs. 6 and 7). 

The Prandtl number functions of figure 10 are similar to those 
predicted by Kaeppeler and Krause (ref. 30) . Their calculations were 
based on the low value for nitrogen dissociation energy and so the pres- 
ent estimates exhibit more distinctly separate maxima. This is because 
of the separation between the dissociation reactions which occurs when 
the higher value for nitrogen dissociation energy is used. Greifinger 
(ref. 31) has also made some estimates for the transport properties of 
air at high temperatures. He assumes constant collision cross sections 
in a manner such that, up to the point where ionization begins, his calcu- 
lated Prandtl number is a constant lacking the structure shown in fig- 
ure 10. However, reference 31 predicts the same order of effect due to 
ionization which has been discussed here. 
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Recall that the nitric oxide in air has teen neglected and that the 
calculated thermodynamic properties of air are not greatly affected by 
this omission. It is not obvious that the transport properties of air 
will be similarly independent of NO formation, since NO is ionized 
rather easily compared to N 2 and 0 2 molecules. Practically all of the 
electrons which appear in air at intermediate temperatures come from 
the NO ionization reaction. However, Gilmore finds that the quantity 
of such electrons is very small, the order of 0.01 percent or less (ref. 9) > 
and from equations (94) and (95) it can be judged that this amount is not 
yet sufficient to alter greatly the transport properties of air. 

At present there are no direct measurement s of the transport 
properties for air at the high temperatures considered here. One experi- 
mental comparison which can be made at the moment is to insert the calcu- 
lated coefficients into a theory for the heat transfer to a simple 
aerodynamic shape, and then compare the results with measured heat- 
transfer rates. For example, reference 32 compares the values for heat 
flux to the stagnation region of a blunt body in high-energy air flow 
which were measured by Rose and Riddell (ref. 33) with the theoretical 
predictions based on the transport coefficients developed in the present 
report. It is found that the calculated heat flux deviates less than 
10 percent from the mean of the measured data. This is a favorable result 
in view of the approximations involved in the theory and of the scatter 
in the experimental data (the order of +20 percent) . It should be empha- 
sized that this comparison is not a sufficient test for the reliability 
of the calculated transport coefficients, as it is possible that errors 
in the coefficients compensate for each other or for approximations in 
the heat-transfer analysis. Nevertheless, the favorable comparison is 
reassuring inasmuch as it means that at least a necessary condition for 
the correctness of the theories is fulfilled. 


CONCLUDING REMARKS 


In conclusion, the thermodynamic properties (including compressibility, 
energy, entropy, specific heat, and the speed of sound) and the transport 
properties (the coefficients of viscosity and thermal conductivity and 
the Prandtl number) have been estimated for air at high temperatures. 

These estimates were made from approximations in closed form, and it is 
found that these approximations give fairly accurately the thermodynamic 
properties over the range of pressures and temperatures for flight through 
the atmosphere. The transport properties were calculated by methods which 
have given reasonably good results for gases at ordinary temperatures and 
for gases which dissociate and react chemically. Moreover, the values 
calculated for the transport properties can be used to predict the same 
order of heating rates that have been measured at the stagnation region 
of bl un t bodies in high-energy air flow. Therefore it is concluded that 
the results presented in this paper are useful engineering approximations 
for the properties of high-temperature air. In particular, it is an 
advantage that the solutions appear in closed form because the functional 
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relationships can be visualized. In addition the air properties can 
conveniently be computed in very small intervals to give tables that can 
be used in obtaining solutions to real air flow by the method of char- 
acteristics, or for numerical solutions of boundary-layer flow. It is 
expected that the calculated coefficients of viscosity and thermal con- 
ductivity will prove to be lower bounds to more precise calculations or 
measurements. The prediction that the Prandtl number for highly ionized 
flow will be small compared to unity is of interest, for it implies that 
boundary- layer regions in such a gas will be very transparent to heat 
flux. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Nov. 18, 1957 
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TABLE I.- PARTITION FUNCTION CONSTANTS FOR THE MAJOR COMPONENTS OF AIR 
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TABLE II.- THERMODYNAMIC PROPERTIES FOR THE MAJOR COMPONENTS OF AIR 


T, 

°K 

(a) Logarithm of the partition function 
for a standard state of one 
atmosphere, Zn Q p 

E— E 

(b) Dimensionless energy, 0 

RT 

(c) Dimensionless specific 

heat, i — 

9 R dT 


Particle 

Particle 

Particle 

n 2 

0 2 

N 

0 

N+ 

0 + 

e“ 

N 2 

°2 

N 

0 

N + 

0 + 

e" 

n 2 

0 2 

N 

0 


0 + 

e“ 

500 

1,000 

1,500 

2,000 

2,500 

3.000 

3.500 

4.000 

4.500 

5.000 

5.500 

6.000 

6.500 
7,000 

7.500 
8,000 

8.500 
9,000 

9.500 
10,000 

10.500 
11,000 

11.500 
12,000 

12.500 

13.000 

13.500 

14.000 

14.500 

15.000 

21.33 

23.79 

25.28 

26.39 

27.26 

27.99 

28.62 

29.17 

29.66 

30.10 
30.50 
30.87 
31.21 
31.52 
31.82 

32.10 
32.36 
32.60 
32.84 
33.06 

22.97 

25.49 

27.05 

28.20 

29.12 
29.88 

30.54 

31.12 
31.64 
32.10 

32.54 
32.93 
33.30 
33-64 

20.68 

21.24 

21.70 

22.08 

22.42 
22.72 
22.98 

23.23 
23.46 

23.67 

23.86 
24.05 

24.23 
24.40 
24.57 
24.72 

24.87 
25.02 
25.16 
25.29 

25.42 
25.55 

25.67 
25.78 
25.90 
26.01 
26.12 

19.86 

20.90 

21.64 
22.21 

22.67 
23.06 
23.40 
23.70 
23.97 
24.21 
24.43 

24.64 

24.83 

25.01 
25.18 

25.34 
25.49 
25.63 

25.76 

25.89 

26.01 
26.13 
26.25 

26.35 
26.46 
26.56 

26.68 
26.75 

26.84 

23.77 

24.01 

24.23 

24.44 

24.63 

24.81 
24.98 
25.14 
25.29 
25.43 
25.57 
25.70 

25.82 
25.94 
26.05 
26.16 
26.27 
26.37 
26.46 
26.56 
26.65 

23.18 

23.41 

23.63 

23.84 

24.02 
24.20 
24.37 
24.53 
24.68 
24.82 
24.96 
25.10 
25.22 
25.35 
25.47 
25.58 
25.70 

25.81 

25.91 

26.02 
26.12 

7.06 
7.30 

7.52 

7.72 
7.90 

8.07 

8.24 
8.39 

8.53 
8.66 
8.79 
8.92 
9.03 
9.14 

9.25 
9.35 
9.45 

9.54 
9.63 

9.72 
9.81 

2.51 

2.62 

2.76 

2.88 

2.97 

3.04 

3.09 

3-14 

3.17 

3.20 

3.22 

3.24 
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3-28 
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3.31 
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3.70 
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1.50 
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1.53 
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2.02 

2.04 

2.07 

2.09 

2.11 

2.13 

2.14 
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2.16 

1.60 
1.57 
1.55 
1.54 
1.54 
1.54 
1.54 

1.54 

1.55 

1.56 

1.57 

1.58 
1-59 
1.60 
1.60 

1.61 
1.62 

1.63 

1.64 

1.64 

1.65 

1.66 
1.66 
1.67 

1.67 

1.68 
1.68 
1.68 
1.68 

1.56 

1.56 

1.57 

1.58 

1.59 

1.60 
1.61 
1.62 

1.63 

1.64 

1.65 

1.66 
1.66 

1.67 

1.68 
1.68 

1.69 

1.70 

1.70 
1-71 

1.71 

1.51 

1.52 

1.53 

1.54 
1.56 
1.58 
1.60 
1.63 
1.66 
1.69 
1.72 
1.75 
1.78 
1.81 
1.84 
1.87 
1.90 
1-93 
1.95 
1.98 
2.00 

1. 

\ 

50 

' 

2.55 

2.92 

3.16 

3.29 

3.36 

3.40 

3.42 

3.44 

3.45 

3.46 
3-47 

3.47 

3.48 
3-48 
3.48 

3.48 

3.49 
3.49 
3.49 
3-49 

2.72 

3.16 

3.37 

3.46 

3.62 

3.71 

3.80 

3.88 

3.94 

3.97 

4.00 

4.00 

3.99 

3-98 

3.96 

3.96 

3.94 

3.92 

3.90 

3.88 

1.50 

1.51 

1.52 

1.56 
1.62 
1.71 
1.82 
1.94 
2.07 
2.19 
2.30 
2.40 
2.49 
2.55 
2.61 

2.64 
2.66 
2.67 
2.67 

2.65 
2.64 
2.61 

2.57 
2.55 

2.52 
2.48 
2.45 

1.52 

1.51 

1.50 

1.51 

1.52 
1.54 
1.56 
1.59 
1.62 
1.65 
1.68 
1.70 
1.72 
1.74 

1.76 

1.77 

1.77 

1.78 
1.78 
1.78 
1.78 
1.78 
1.78 
1.78 
1.78 
1.77 
1.77 
1.76 

1.76 

1.63 

1.66 

1.69 

1.72 

1.74 

1.75 

1.77 

1.78 
1-79 

1.80 

1.81 
1.82 
1.82 
1.83 

1.83 

1.84 
1.84 

1.84 

1.85 
1.85 
1.85 

1.57 

1.61 

1.67 

1.74 

1.82 

1.91 

2.00 

2.09 

2.18 

2.26 

2.34 

2.41 

2.46 

2.52 

2.56 

2.60 

2.63 

2.65 

2.66 
2.66 
2.66 

1.! 

\ 

>0 

f 
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TABLE III.- LOGARITHMS OF THE CHEMICAL EQUILIBRIUM CONSTANTS AND THEIR DERIVATIVES 

FOR THE MAJOR REACTIONS IN AIR 


4 =- 

o\ 


T, 

°K 

(a) to Kp 

. . d to Kt, 

(b) T p 

dT 

Reaction 

Reaction 

0 2 20 

N 2 2R 

0 0 + +e~ 

N I^+e" 

0 2 20 

N 2 -> 2N 

0 0 + +e- 

N I^+e- 

500 

1,000 

1.500 
2,000 

2.500 

3.000 

3.500 

4.000 

4.500 

5.000 

5.500 

6.000 

6.500 

7.000 

7.500 

8.000 

8.500 
9,000 

9.500 
10,000 

10.500 
11,000 

11.500 
12,000 

12.500 

13.000 

13.500 

14.000 

14.500 

15.000 

-44.78 
-24.58 
-14.42 
-8.30 
-4.20 
-1.27 
• 94 
2.65 
4.03 
5.16 
6.10 
6.90 
7-59 
8.19 
8.72 
9-18 
9.60 
9*98 
IO.32 

-41.62 

-30.06 

-22.33 

-16.80 

-12.63 

-9-38 

-6.77 

-4.62 

-2.82 

-1.29 

.03 

1.19 

2.22 

3.13 

3.95 

4.69 

5.37 

-25.33 
-22.25 
-19.62 
-17.40 
-15-48 
-13.80 
-12.32 
-11.01 
-9.83 
-8.77 
-7- 81 
-6.93 
-6.12 
-5-38 
-4.70 
-4.06 
-3-47 
-2.91 
-2.40 
-1.91 
-1.45 

-25.92 

-22.61 

-19.84 

-17-48 

-15.45 

-13.68 

-12.12 

-10.74 

-9.50 

-8.39 

-7.39 

-6.48 

-5.65 

-4.89 

-4.19 

-3.54 

-2.94 

-2.38 

-1.86 

-1.37 

-.91 

60.43 

40.54 

30.56 

24.53 

20.51 

17.62 

15.45 

13.76 

12.41 

11.31 

10.40 

9-63 

8.97 

8.40 

7.91 

7-48 

7.10 

6.76 

6.46 

57.72 

46.31 
38.70 
33.27 
29.20 
26.05 
23.55 
21.53 

19.87 
18.47 

17.31 
16.30 
15.43 
14.67 
14.00 
13.41 

12.88 

34.06 

31.19 
28.79 
26.77 
25.04 
23-55 
22.25 
21.10 
20.09 

19.19 

18.38 
17.65 
16.99 

16.39 
15.84 

15.34 
14.88 
14.46 

14.06 
13.69 

13.35 

36.26 

33.17 

30.59 
28.39 
26.51 

24.86 

23.42 

22.15 
21.01 
20.00 
19. 08 

18.25 

17.49 

16.81 

16.18 

15.60 
15.07 
14.57 
14.12 
13.69 
13.30 
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TABLE IV.- THERMODYNAMIC PROPERTIES FOR AIR 


°K 


(a) Compressibility, 

z 


(b) Dimensionless energy, ~ 


(c) Dimensionless entropy, 

ZS 

R 

Pressure, atmospheres 

Pressure, atmospheres 

Pressure, atmospheres 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 


500 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

2.52 

2.52 

2.52 

2.52 

2.52 

2.52 

2.52 

21.1 

23.9 

25.7 

28.0 

30.3 

32.6 

34.9 


§ 1,000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

2.65 

2.65 

2.65 

2.65 

2.65 

2.65 

2.65 

23.7 

26.0 

28.3 

30.6 

32.9 

35.2 

37.5 


2 1,500 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

2.80 

2.80 

2.80 

2.80 

2.80 

2.80 

2.80 

25.3 

27.6 

29.9 

32.2 

34.5 

36.8 

39.2 


5 2,000 

1.000 

1.000 

1.000 

1.001 

1.002 

1.005 

1.016 

2.92 

2.92 

2.92 

2.93 

2.96 

3.06 

3.40 

2 6.6 

28.9 

31.2 

33.5 

35.8 

38.2 

40.9 

1 

§ 2,500 

1.000 

1.001 

1.003 

1.011 

1.033 

1.087 

1.163 

3.02 

3.03 

3.09 

3.27 

3.78 

5.05 

6.84 

27.5 

29.9 

32.2 

34.7 

37.6 

41.4 

45.8 


a 3,000 

1.003 

1.009 

1.026 

1.072 

1.149 

1.192 

1.201 

3.13 

3.25 

3.58 

4.48 

5.98 

6.82 

7.00 

28.4 

30.8 

33.5 

36.9 

41.0 

44.6 

1+7.6 

Sd 3,500 

1.012 

1.035 

1.091 

1.152 

1.196 

1.204 

1.211 

3.33 

3.71 

4.65 

5.65 

6.39 

6.57 

6.81 

29.3 

32.0 

35.4 

39.1 

42.6 

45.6 

48.6 

o 

O ^.000 

1.033 

1.088 

1.164 

1.195 

1.209 

1.228 

1.287 

3.66 

4.46 

5.56 

6.00 

6.33 

6.86 

8.52 

30.2 

33.5 

37.3 

40.5 

43.6 

46.9 

51.5 


g) ^,500 

1.070 

1.147 

1.191 

1.214 

1.244 

1.336 

1.575 

4.12 

5.10 

5.66 

6.13 

6.88 

9.19 

15.17 

31.3 

34.9 

38.2 

ki.5 

45.1 

50.4 

59.1* 


5 5,000 

1.115 

1.180 

1.217 

1.252 

1.358 

1.620 

1.908 

4.58 

5.32 

5.91 

6.71 

9.10 

15.00 

21.5 

32.3 


39.2 

42.9 

48.4 

58.0 

68.8 


° 5,500 

1.153 

1.215 

1.248 

1.347 

1.597 

1.896 

1.987 

4.86 

5.66 

6.33 

8.37 

13.51 

19.63 

21.5 

33.3 

36TF 

40.8 

45.5 

54.2 

64.7 

71.2 


6,000 

1.176 

1.237 

1.316 

1.527 

1.847 

1.97? 

2.010 

4.96 

5.89 

7-37 

11.36 

17.39 

19.90 

20.7 

33.9 

37.7 

42.2 

49.6 

59.9 

66.9 

72.2 


6,500 

1.226 

1.280 

1.436 

1.750 

1.958 

2.010 

2.032 

5.51 

6.45 

9.18 

14.66 

18.30 

19.31 

19.9 

35.6 

38.9 

44.9 

54.3 

62.5 

68.1 

73.4 


§ 7,000 

1.250 

1.351 

1.605 

1.900 

1.988 

2.028 

2.088 

5.76 

7.40 

11.30 

16.36 

17' -7? 

18.69 

20.2 

35.8 

40.5 

48.3 

57.4 

442 

69.1 

75.3 


S 7,500 

1.287 

1.456 

1.775 

1.940 

2.022 

2.068 

2.210 

6.18 

8.75 

13.63 

16.16 

17.59 

18.67 

22.0 

36.8 

42.7 

51.6 

58.6 

64.7 

70.5 

78.7 


* 8,000 

1.343 

1.589 

1.891 

1.986 

2.046 

2.144 

2.446 

6.82 

10.37 

14.73 

16.11 

17.32 

19.49 

26.2 

38.0 

45.I 

53.8 

59-8 

65.7 

72.5 

84.6 


§ 8,500 

1.417 

1.722 

1.950 

2.028 

2.090 

2.284 

2.827 

7-65 

11.83 

14.94 

16.21 

17.50 

21.6 

32.9 

39-1 

47.4 

55-1 

SlTcT 

67.O 

76.1 

93.2 

H 

i— i 

S 9,000 

1.510 

1.832 

1.990 

2.052 

2.166 

2.510 

3.282 

8.67 

12.86 

14.95 

16.08 

18.36 

25.2 

40.4 

4i.o 

49.4 

56.0 

61.9 

69.0 

81.1 

103.0 

9,500 

1.612 

1.904 

2.020 

2.092 

2.286 

2.834 

3.646 

9.67 

13.30 

15.02 

16.22 

19.89 

30.2 

45.5 

42.7 

50.7 

57.0 

63.0 

71.7 

87.8 

110.6 

w 

aJ 

a 10,000 

1.713 

1.970 

2.048 

2.150 

2.462 

3.202 

3.844 

10.60 

1?.6? 

14.86 

16.71 

22.3 

35.6 

47.1 

44.4 

51.? 

57.8 

64.5 

75.3 

95.0 

114.8 

o 

£,10,500 

1.799 

2.024 

2.074 

2.236 

2.700 

3-528 

3.932 

11.29 

13.75 

14.86 

17.61 

25.3 

39.8 

46.7 

46.4 

53.1 

58.6 

66.4 

79.9 

100.9 

116.8 


gii,ooo 

1.864 

2.036 

2.112 

2.352 

2.984 

3-746 

3.969 

11.70 

13-78 

15.04 

18.97 

29.0 

41.8 

45.5 

47.5 

53.6 

59.6 

68.6 

85.0 

105.4 

118.0 


^11, 500 

1.910 

2.052 

2.176 

2.506 

3.274 

3.868 

3.986 

11.89 

13.66 

15.50 

20.8 

32.9 

42.2 

44.0 

^7.7 

54.2 

61.1 

71.4 

90.1 

107.8 

118.8 


» 12,000 

1.941 

2.074 

2.234 

2.696 

3.522 

3.930 

3.993 

11.95 

13.63 

16.05 

23.O 

35.5 

41.7 

42.6 

48.3 

54.8 

62.1 

74.7 

94.4 

109.2 

119.3 


12,500 

l.?96 

2.102 

2.320 

2.912 

3.700 

3.962 

3.996 

12.16 

13.68 

16.86 

25.5 

36.9 

40.8 

41.2 


55.5 

63.7 

78.3 

97.4 

110.1 

119.8 

1— 1 

§13,000 

2.021 

2.138 

2.426 

3.136 

3. 818 

3-980 

3.998 

12.30 

13.85 

17.90 

27.9 

37.4 

39.7 

39.9 

W^9 

56.3 

65.5 

81.8 

99.5 

110.8 

120.2 

H 
1 — ( 

^13,500 

2.058 

2.182 

2.554 

3.348 

3.890 

3-988 

3.999 

12.44 

14.12 

19.16 

29.9 

37.3 

38.6 

38.7 

50.6 

57.2 

67.6 

85-1 

100.9 

111.3 

120.7 

<D 

g 14,000 

2.074 

2.234 

2.700 

3.526 

3.932 

3-994 

3.999 

12.49 

14.48 

20.6 

31.4 

36.7 

37.5 

37.6 

51.2 

58.1 

69.8 

87-8 

101.8 

111.8 

121.1 

co 

a] 

311,500 

2.102 

2.320 

2.912 

3.700 

3.962 

3.996 

4.000 

12.48 

15.24 

22.7 

32.8 

36.1 

36.5 

36.6 

51.7 

59.5 

72.8 

90.4 

102.6 

112.2 

121.5 

O 

015,000 

2.120 

2.374 

3.030 

3.770 

3.998 

3-998 

4.000 

12.44 

15.56 

23.6 

32.8 

35.6 

35.6 

35.6 

52.2 

60.4 

74.6 

91.6 

103.4 

112.6 

121.9 


TABLE TV.- THERMODYNAMIC PROPERTIES FOR AIR - Concluded 


£ 


T, 

°K 

(d) Dimensionless specific heat at 
constant density, i ) 

(e) Dimensionless specific heat at 

constant pressure, — ) 

R \dT/p 

(f) Dimensionless speed of sound 
parameter, ^72 

Pressure, atmospheres 

Pres sure , atmo sphe re s 

Pressure, atmospheres 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 


500 

2.59 

2.59 

2.59 

2.59 

2.59 

2.59 

2.59 

3.59 

3.59 

3.59 

3.59 

3.59 

3.59 

3.59 

1.39 

1.39 

1.39 

1.39 

1.39 

1.39 

1.39 

d 

1,000 

2.96 

2.96 

2.96 

2.96 

2.96 

2.96 

2.96 

3.96 

3.96 

3.96 

3.96 

3.96 

3.96 

3.96 

1.34 

1.34 

1.34 

1.34 

1.34 

1.34 

1.34 

•H 

1,500 

3.20 

3.20 

3.20 

3.20 

3.20 

3.20 

3.20 

4.20 

4.20 

4.20 

4.20 

4.20 

4.20 

4.20 

1.31 

1.31 

1.31 

1.31 

1.31 

1.31 

1.31 


2,000 

3.33 

3.33 

3.33 

3.54 

3.98 

5.49 

10.12 

4.33 

4.33 

4.33 

4.56 

5.03 

6.65 

11.66 

1.30 

1.30 

1.30 

1.29 

1.25 

1.21 

1.14 

•H 

a 

2,500 

3.55 

3.69 

4.37 

6.40 

n.60 

20.7 

17.45 

4.56 

4.71 

5-37 

7-68 

13.47 

23.9 

20.2 

1.29 

1.28 

1.23 

1.19 

1.14 

1.12 

1.13 

1— 1 O 

CO 

3,000 

3.94 

5.05 

7.99 

14.10 

14.97 

6.25 

3.75 

5.01 

6.24 

9-55 

16.60 

17.68 

7.77 


1.27 

1.23 

1.1 8 

1.15 

1.15 

1.24 

1.53 

<D CO 
tQ -H 

3,500 

5.05 

7.82 

12.32 

11.55 

4.61 

5.23 

8.99 

6.27 

9.45 

14.79 

13.94 

5.97 

6.^5 

10.58 

1.23 

1.19 

1.17 

1.18 

1.29 

i.25 

1.17 

cd 

0 

4,000 

6.64 

10.06 

8.64 

4.33 

6.96 

14.16 

35.9 

8.18 

12.32 

10.68 

5.67 

8.44 

16.30 

40.7 

1.21 

1.19 

1.21 

1.30 

1.2l 

1.14 

1.10 

d 

a) 

4,500 

8.00 

8.45 

4.66 

7.74 

16.87 

42.2 

85.1 

9.93 

10.60 

6.06 

9.34 

19.46 

48.7 

101.2 

1.21 

1.22 

1.29 

17201 

1.13 

1.10 

1.09 

| 

5,000 

7.94 

5.46 

7.56 

16.20 

39.3 

72.3 

40.6 

10.04 

7.04 

9. 18 

18. 90 

46.0 

87-0 

47.4 

1.22 

1.27 

1.21 

1.14 

l.H 

1.10 

1.12 

0 

5,500 

6.57 

6.61 

13.27 

31.4 

59.7 

37-5 

9.39 

8.46 

8.16 

15.69 

37.1 

72.4 

44.4 

11.97 

1.24 

1.23 

1.16 

1.12 

1.11 

1.13 

1.27 


6,000 

5.11 

9.91 

22.5 

46.9 

40.8 

10.60 

8.27 

6.68 

11.96 

26.8 

57.2 

49.2 

13.37 

10.59 

1.29 

1.19 

1.14 

1.12 

1.14 

1.24 

1.28 


6,500 

7.36 

15.03 

33.5 

44.0 

15.96 

7.86 

16.02 

9- 08 

18.04 

40.8 

54.0 

19.57 

IO.17 

19.07 

1.22 

1.17 

1.14 

1.14 

1.20 

1.29 

1.18 

d 

7,000 

9.96 

21.6 

39.5 

24.7 

7.44 

13.36 

31.7 

12.17 

26.2 

49.2 

30.2 

9.84 

16.18 

36.6 

1.20 

1.16 

1.14 

1.17 

1.31 

1.20 

1.13 

0 

•H 

7,500 

13.17 

28.0 

33.8 

16.57 

10.45 

23.O 

60.1 

16.09 

34.8 

42.1 

20.5 

13.01 

27.O 

69.8 

1.18 

1.15 

1.16 

1.20 

TT2JT 

1.15 

1.12 

Is 

8,000 

17.H 

31.5 

21.8 

7.46 

15.92 

38.8 

100.6 

21.2 

40.0 

27.2 

9.88 

19.16 

45.4 

121.2 

1.18 

1.16 

1.19 

1.32 

1.19 

1.13 

1.11 

•H 

O 

8,500 

21.0 

29.5 

13.08 

n.05 

24.1 

62.5 

139.4 

26.4 

37-5 

16.52 

13.73 

28.5 

74.4 

173.3 

1.18 

1.17 

1.23 

T7l3 

1.16 

1.12 

1.11 

O 

H CO 

9,000 

23.9 

23.1 

6.52 

15.06 

36.3 

90.7 

139.0 

30.6 

29.4 

8-79 

18.29 

43.0 

111.0 

172.8 

1.18 

1.19 

1.34 

1.21 

1.14 

1.12 

1.12 

M CO 
•H 

9,500 

25.0 

16.34 

8.52 

20.9 

52.1 

114.5 

94.1 

32.2 

20.7 

10.96 

25;2 

62.8 

144.1 

114.6 

1.18 

1.22 

TTSF 

1.18 

1.14 

1.12 

1.13 

<d 

CO 

10,000 

23.4 

8.77 

12.49 

28.3 

70.0 

116.2 

49.0 

30.6 

11.51 

15.49 

34.0 

86.3 

146.7 

59.4 

1.20 

1.29 

1.23 

1.16 

1.13 

1.13 

1.17 

a3 d 
O <D 

10,500 

20.3 

8.54 

15.79 

37.5 

86.5 

9 2.0 

25.1 

26.2 

11.02 

19.40 

45.6 

109.5 

114.8 

31.2 

1.21 

1725" 

1.21 

1.16 

1.13 

1.14 

1.24 

bO 

O 

11,000 

16.33 

9.76 

19.82 

48.1 

95.1 

59.6 

15.14 

21.2 

12.44 

24.2 

59.3 

121.7 

73-3 

20.2 

1.23 

1.26 

1.19 

1.15 

1.14 

1.17 

1.32 

-P 

11,500 

12.80 

11.40 

24.8 

59.1 

90.3 

35.4 

10.81 

16.71 

14.35 

30.2 

74.0 

115.1 

43.8 

15.30 

1.25 

1.24 

1.18 

1.15 

1.15 

1.20 

1.41 

•H 

12,000 

10.15 

13.37 

30.3 

68.4 

73-6 

21.9 

8.76 

13.36 

16.74 

37.4 

87.5 

93.1 

28.2 

13.00 

1.28 

1.23 

1.17 

1.15 

1.16 

1.26 

1.48 


12,500 

5.27 

15.47 

36.1 

73.6 

53.4 

14.20 

7-82 

7.34 

19.35 

45.1 

95.8 

67.2 

19.26 

11.93 

1.37 

1.22 

1.17 

1.16 

1.18 

1.34 

1.52 

H g 

13,000 

^82 

18.06 

42.3 

73.5 

36.8 

10.82 

7.39 

9.15 

22.7 

53-9 

96.3 

46.6 

15.40 

11.44 

1.34 

1.21 

1.17 

1.17 

1.21 

1.42 

1.55 

H *H 
M -P 

13,500 

10.05 

20.9 

47.9 

66.8 

24.8 

9.16 

7.01 

13.03 

26.3 

62.2 

87.5 

32.0 

13.50 

11.00 

1.28 

1.21 

1.18 

1.13 

1.26 

1.47 

1.57 

ccJ 

to rsJ 

14,000 

11.14 

24.0 

53-1 

56.8 

17.73 

8.38 

7.02 

14.33 

30.3 

69.5 

73-5 

23.6 

12.59 

11.02 

1.26 

1.20 

1.18 

1.19 

1.31 

1.50 

1.57 

2-d 14,500 

12.81 

28.6 

57.5 

42.3 

U.95 

7.66 

7.02 

16.41 

36.5 

76.0 

54.4 

17.42 

H.76 

11.02 

1.25 

1.20 

1.18 

1.21 

1.43 

1.53 

1.57 

0 0 

M 

15,000 

13.55 

30.5 

56.8 

34.4 

7.63 

7-33 

7.03 

17.50 

39.5 

76.1 

44.6 

11.73 

11.38 

11.03 

1.26 

1.20 

1.19 

1.24 

1.54 

1.55 

1.57 
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TABLE V.- COLLISION CROSS SECTIONS 


T, 

°K 

So, in 

10" 16 cm2 

s(n 2 -n) 

S(N-N) 

S(N-e) 

S(e-e) 

S’ (n 2 -n) 

S' (N-N) 

So 

So 

So 

So 

So 

s o 

500 

38.4 

0.946 

0.894 



0.877 

0.761 

1,000 

34.9 

.920 

.838 



.843 

.703 

1,500 

33.7 

.889 

• 785 



.817 

.652 

2,000 

33.2 

.886 

.742 



• 79^ 

.611 

2,500 

32.8 

.846 

.705 



• 775 

.578 

3,000 

32.6 

.830 

.676 



.759 

• 551 

3,500 

32.4 

.815 

.650 



.745 

.527 

4,000 

32.3 

.803 

.628 



• 733 

.507 

4,500 

32.2 

.792 

.608 



.722 

.489 

5,000 

32.1 

.782 

.591 



.712 

.473 

5,500 

32.0 

• 773 

.575 

0.397 

89.9 

.703 

.458 

6,000 

32.0 

.764 

.561 

• 380 

75.6 

.695 

.445 

6,500 

31.9 

• 757 

.548 

.366 

64.5 

.688 

.433 

7,000 

31.9 

.750 

.536 

.353 

55-7 

.681 

.422 

7,500 

31.9 

• 743 

.524 

.342 

48.6 

.674 

.412 

8,000 

31.8 

• 737 

.514 

.331 

42.8 

.66 8 

.402 

8,500 

31.8 

• 731 

.504 

.321 

37.9 

.662 

• 393 

9,000 

31.8 

.725 

.495 

.313 

33.8 

.657 

• 385 

9,500 

31.8 

.720 

.486 

.304 

30.4 

.652 

• 377 

10,000 

31.8 

.715 

.478 

.297 

27.4 

.647 

• 370 

10,500 

31.7 

.710 

.470 

.290 

24.9 

.642 

.363 

11,000 

31.7 

.706 

.463 

.283 

22.7 

.637 

.356 

11,500 

31.7 

.701 

.456 

.281 

20.8 

.633 

.350 

12,000 

31.7 

.697 

.448 

.270 

19.09 

.629 

.342 

12,500 

31.7 

.693 

.443 

.266 

17.60 

.625 

.338 

13,000 

31.7 

.689 

.437 

.261 

16.27 

.621 

• 332 

13,500 

31.7 


.431 

.256 

15.10 


.327 

14,000 

31.7 


.426 

.252 

14.04 


.322 

14,500 

31.6 


.420 

.247 

13.09 


.316 

15,000 

31.6 


.415 

.243 

12.24 


.312 


S(N 2 -N) = S(0 2 -0) = s(n 2 -o) 

S(N-N) = S(O-O) = S(N-O) = StN-N*) = S(0-0 + ) 
S(N-e-) = S(O-e-) 

S' (N 2 -N) = S' (0 2 -0) = S* (N 2 -0) 

S*(N-N) = S' (O-O) = S' (N-0) = S' (N-N + ) = S' (0-0 + ) 
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TABLE VI.- TRANSPORT PROPERTIES OF AIR 


(a) Coefficient of viscosity 


T, 

°K 



Ratio 

nAio 



Reference 
coefficient, 
tfo, in 1° ' 6 

Pre s sure , atmo sphere s 

lb sec 

gm 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

ft 2 

cm-sec 

500 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

0.558 

267 

1,000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

.868 

416 

1,500 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.100 

527 

2,000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.293 

619 

2,500 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.461 

700 

3,000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.612 

772 

3,500 

1.000 

1.001 

1.003 

1.006 

1.010 

1.010 

1.011 

1.751 

838 

4,000 

1.003 

1.008 

1.016 

1.020 

1.022 

1.024 

1.032 

1.879 

899 

4,500 

1.010 

1.022 

1.029 

1.033 

1.038 

1.055 

1.096 

1.999 

957 

5,000 

1.022 

1.036 

1.043 

1.051 

1.074 

1.128 

1.181 

2.11 

1011 

5,500 

1.036 

1.052 

1.060 

1.O86 

1.146 

1.209 

1.227 

2.22 

1062 

6,000 

1.050 

1.067 

1.090 

1.148 

1.228 

1.257 

1.259 

2.32 

1112 

6,500 

1.072 

1.090 

1.139 

1.229 

1.276 

1.289 

1.277 

2.42 

1159 

7,000 

1.089 

1.124 

1.208 

1.294 

1.317 

1.309 

1.273 

2.52 

1204 

7,500 

1.112 

1.175 

1.283 

1.332 

1.342 

1.315 

1.221 

2.61 

1247 

8,000 

1.143 

1.238 

1.342 

1.371 

1.355 

1.291 

1.086 

2.69 

1289 

8,500 

1.185 

1.307 

1.386 

1.392 

1.354 

1.220 

.841 

2.78 

1330 

9,000 

1.238 

1.368 

1.425 

1.405 

1.328 

1.085 

.534 

2.86 

1370 

9,500 

1.298 

1.418 

1.447 

1.404 

1.267 

.871 

.280 

2.94 

1408 

10,000 

1.361 

1.468 

1.460 

1.389 

1.162 

.615 

.1397 

3.02 

1446 

10,500 

1.418 

1.497 

1.467 

1.352 

1.004 

.384 

.0788 

3.10 

1482 

11,000 

1.467 

1.515 

1.464 

1.287 

.806 

.226 

.0556 

3.17 

1518 

11,500 

1.509 

1.529 

1.450 

1.190 

• 597 

.1362 

.0466 

3.24 

1552 

12,000 

1.549 

1.543 

1.425 

1.065 

.415 

.0930 

.0435 

3.31 

1586 

12,500 

1.577 

1.542 

1.376 

.908 

.279 

.0688 

.0436 

3.38 

1620 

13,000 

1.592 

1.539 

1.312 

.741 

.1918 

.0598 

.0453 

3-45 

1652 

13,500 

1.617 

1.528 

1.230 

.580 

.1363 

.0569 

.0471 

3.52 

1684 

14,000 

1.628 

1.508 

1.129 

.443 

.1061 

.0572 

.0506 

3-58 

1716 

14,500 

1.631 

1.460 

• 968 

.306 

.0843 

.0576 

.0543 

3-65 

1747 

15,000 

i.64o 

1.436 

.882 

.254 

.0615 

.0598 

.0581 

3.71 

1777 
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TABLE VI.- TRANSPORT PROPERTIES OF AIR - Continued 


(b) Coefficient of thermal conductivity 


T, 

°K 

Ratio k/k Q 

Reference 
coefficient, 
k Q , in 10 -6 

Pre s sure , atmo sphere s 

Btu 

watt 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

ft sec °R 

cm °K 

500 

1.021 

1.021 

1.021 

1.021 

1.021 

1.021 

1.021 

5.84 

364 

1,000 

1.100 

1.100 

1.100 

1.100 

1.100 

1.100 

1.100 

9.10 

567 

1,500 

1.150 

1.150 

1.150 

1.150 

1.150 

1.150 

1.150 

11.53 

719 

2,000 

1.177 

1.177 

1.177 

1.251 

1.460 

2.09 

3-99 

13.55 

844 

2,500 

1.256 

1.317 

1.619 

2.50 

4.63 

7.67 

5.50 

15.31 

954 

3,000 

1.421 

1.928 

3-20 

5.48 

5.02 

2.19 

1.465 

16.90 

1053 

3,500 

1.941 

3.15 

4.72 

3.96 

1.719 

2.11 

3.71 

18.35 

1143 

4,000 

2.69 

3.94 

2.99 

1.600 

2.91 

6.04 

15.03 

19.69 

1227 

4,500 

3.22 

3-06 

1.714 

3.32 

7.34 

17.65 

30.5 

21.0 

1305 

5,000 

3.07 

1-997 

3.29 

7-18 

16.63 

25.8 

11.84 

22.1 

1379 

5,500 

2.46 

2.91 

5.99 

13.71 

22.2 

11.40 

3.54 

23.3 

1449 

6,000 

1.930 

4.53 

10.19 

18.74 

13.09 

3.96 

7.30 

24.3 

1516 

6,500 

3-35 

6.98 

14.50 

15.39 

5.49 

7.18 

14.69 

25.4 

1580 

7,000 

4.69 

9-97 

15.69 

8.32 

3-28 

12.82 

29.0 

26.4 

1642 

7,500 

6.31 

12.48 

12.24 

5.92 

10.47 

22.3 

54.3 

27.3 

1701 

8,000 

8.21 

13.19 

7-80 

3.42 

16.37 

37.5 

87-6 

28.2 

1759 

8,500 

9.86 

11.55 

5.10 

11.99 

25.0 

59.3 

110.9 

29.1 

1814 

9,000 

10.90 

8.79 

3.26 

16.83 

37-5 

82.4 

97-5 

30.0 

1868 

9,500 

10.88 

6.38 

10.09 

23.6 

52.8 

96.5 

59-1 

30.8 

1921 

10,000 

9.87 

4.07 

15.44 

32.0 

69.O 

90.3 

30.2 

31.6 

1972 

10,500 

8.33 

11.18 

19-84 

41.8 

81.2 

64.9 

16.54 

32.4 

2020 

11,000 

6.84 

13.44 

25.1 

52.9 

83-7 

40.0 

11.41 

33-2 

2070 

11,500 

5-59 

16.25 

31.2 

63.2 

74.2 

24.9 

9.54 

34.0 

2120 

12,000 

4.79 

19.56 

38.2 

70.9 

57 -6 

17.18 

8.98 

34.7 

2160 

12,500 

3.34 

23.0 

45.0 

73.2 

41.4 

13.07 

9. 08 

35.5 

2210 

13,000 

ll.ll 

27.O 

52.2 

69.8 

29.6 

11.67 

9.50 

36.2 

2250 

13,500 

17.62 

31.3 

58.2 

61.5 

22.0 

11.37 

9.98 

36.9 

2300 

14,000 

20.0 

35.7 

62.6 

51.1 

17.92 

11.61 

10.73 

37.6 

2340 

14,500 

23.3 

42.1 

65.3 

38.7 

15.24 

11.93 

11.51 

38.2 

2380 

15,000 

25.3 

45-3 

64.2 

33.1 

12.72 

12.52 

12.31 

38.9 

2420 
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TABLE VI.- TRANSPORT PROPERTIES OF AIR - Concluded 


(c) Prandtl number 

T, 

°K 



Pressure , 

, atmospheres 


100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

500 

0.738 

0.738 

0.738 

0.738 

0.738 

0.738 

0.738 

1,000 

.756 

• 756 

• 756 

.756 

• 756 

• 756 

• 756 

1,500 

.767 

.767 

• 767 

.767 

.767 

• 767 

.767 

2,000 

• 773 

• 773 

• 773 

.766 

.724 

.668 

.614 

2,500 

.762 

• 751 

.696 

.645 

.611 

.654 

• 771 

3,000 

.740 

.680 

.627 

.636 

.740 

• 745 

.714 

3,500 

.678 

.631 

.66 0 

.744 

• 737 

.6 58 

.606 

4,000 

.640 

.662 

.762 

• 759 

.619 

.580 

• 587 

4,500 

.654 

• 743 

• 752 

.610 

• 578 

.611 

.764 

5,000 

.702 

.767 

.611 

• 581 

.624 

• 799 

• 993 

5,500 

• 748 

.620 

• 583 

.617 

• 785 

• 989 

.871 

6,000 

.763 

• 592 

.602 

• 736 

• 969 

.891 

.384 

6,500 

.610 

• 592 

• 673 

.906 

• 955 

• 383 

.348 

7,000 

• 593 

.620 

• 796 

• 986 

• 830 

.346 

• 337 

7,500 

• 595 

.688 

• 927 

• 96 9 

• 350 

• 334 

• 330 

8,000 

.620 

.788 

• 983 

.648 

• 332 

.328 

.316 

8,500 

.666 

.891 

• 943 

• 335 

• 324 

.321 

.276 

9,000 

• 730 

.961 

.807 

.321 

• 320 

.307 

• 1987 

9,500 

.806 

.966 

• 330 

• 314 

.316 

• 273 

.1140 

10,000 

.886 

.872 

• 308 

.310 

• 313 

.210 

• 0577 

10,500 

• 937 

.310 

.301 

• 309 

.284 

.1427 

.0312 

11,000 

• 955 

.294 

.296 

• 303 

.246 

.0870 

.0207 

11,500 

• 947 

.284 

• 295 

• 293 

.1945 

.0503 

• 0157 

12,000 

.908 

.277 

• 293 

.276 

.1409 

.0321 

.0132 

12,500 

.728 

.272 

.290 

.250 

.0949 

.0213 

.0120 

13,000 

• 275 

.272 

.284 

.215 

.0634 

.0166 

.0115 

13,500 

.251 

.270 

.276 

• 1733 

.0416 

.0142 

.0109 

14,000 

.245 

.269 

.263 

• 1338 

.0293 

.0130 

.0109 

14,500 

.241 

.265 

• 237 

.0903 

.0202 

.0119 

.0109 

15,000 

.238 

.263 

.220 

.0719 

.0119 

.0114 

.0109 


TABLE VII.- PARTIAL COEFFICIENTS 


T, 

(a) Specific heat. 

) 

(b) Thermal conductivity ratio, — 

^o 

°K 

















Pressure, atmospheres 




Pressure, atmospheres 



100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

500 

3.59 

3.59 

3.59 

3.59 

3.59 

3.59 

3.59 

1.021 

1.021 

1.021 

1.021 

1.021 

1.021 

1.021 

1,000 

3.96 

3.96 

3.96 

3.96 

3.96 

3.96 

3.96 

1.100 

1.100 

1.100 

1.100 

1.100 

1.100 

1.100 

1,500 

4.20 

4.20 

4.20 

4.20 

4.20 

4.20 

4.20 

1.150 

1.150 

1.150 

1.150 

1.150 

1.150 

1.150 

2,000 

4.33 

4.33 

4.33 

4.33 

4.33 

4.33 

4.34 

1.177 

1.177 

1.177 

1.177 

1.177 

1.178 

1.181 

2,500 

4.41 

4.41 

4.41 

4.41 

4.42 

4.45 

4.47 

1.195 

1.195 

1.196 

1.198 

1.206 

1.225 

1.252 

3,000 

4.46 

4.47 

4.47 

4.48 

4.51 

4.52 

4.53 

1.205 

1.208 

1.215 

1.235 

1.268 

1.285 

1.289 

3,500 

4.50 

4.52 

4.52 

4.54 

4.56 

4.55 

4.56 

1.218 

1.230 

1.259 

1.289 

1.310 

1.314 

1.320 

4,000 

4.53 

4.55 

4.57 

4.58 

4.58 

4.60 

4.65 

1.238 

1.270 

1.312 

1.329 

1.340 

1.356 

1.406 

4,500 

4.57 

4.59 

4.60 

4.61 

4.64 

4.74 

4.97 

1.269 

1.319 

1.346 

1.366 

1.396 

1.482 

1.674 

5,000 

4.60 

4.61 

4.64 

4.68 

4.81 

5.12 

5.46 

1.3H 

1.357 

1.389 

1.428 

1.541 

1.777 

1.983 

5,500 

4.62 

4.66 

4.71 

4.85 

5.20 

5.62 

5.76 

1.351 

1.406 

1.448 

1.569 

1.829 

2.08 

2.14 

6,000 

4.64 

4.71 

4.85 

5.20 

5.73 

5.94 

6.00 

1.383 

1.457 

1.566 

1.824 

2.14 

2.25 

3.66 

6,500 

4.72 

4.82 

5-12 

5.71 

6.11 

6.21 

6.25 

1.459 

1.543 

1.767 

2.14 

2.33 

3.95 

4.37 

7,000 

4.78 

5.00 

5-53 

6.16 

6.35 

6.44 

6.56 

1.519 

1.685 

2.05 

2.38 

2.47 

4.73 

4.93 

7,500 

4.89 

5-28 

6.02 

6.40 

6.58 

6.68 

6.97 

1.606 

1.889 

2.33 

2.52 

4.99 

5.35 

5.34 

8,000 

5.04 

5.68 

6.42 

6.66 

6. 78 

6.97 

7.57 

1.733 

2.14 

2.55 

2.66 

5.71 

5-87 

5-57 

8,500 

5.27 

6.06 

6.66 

6.91 

6.97 

7.34 

8.41 

1.898 

2.40 

2.70 

5.91 

6.33 

6.26 

5-58 

9,000 

5.55 

6.42 

6.87 

6.98 

7.20 

7.85 

9-35 

2.10 

2.61 

2.83 

6.65 

6.87 

6.52 

5.46 

9,500 

5.85 

6.66 

6.98 

7.H 

7-49 

8.54 

10.11 

2.32 

2.78 

6.39 

7.34 

7-33 

6.6 3 

5.47 

10,000 

6.15 

6.89 

7.07 

7.26 

7-86 

9.29 

10.52 

2.53 

2.93 

7.52 

7-95 

7.69 

6.66 

5.72 

10,500 

6.41 

7.02 

7.12 

7.45 

8.34 

9-95 

10.73 

2.71 

7.32 

8.28 

8.50 

7.96 

6.76 

6.15 

11,000 

6.59 

7.05 

7.20 

7.67 

8.91 

10.41 

10.84 

2.86 

8.16 

9.00 

9.01 

8.16 

7.02 

6.69 

11,500 

6.71 

7.06 

7-33 

7-97 

9.49 

IO.67 

10.91 

2.97 

9.01 

9.67 

9.45 

8.36 

7.41 

7-28 

12,000 

6.76 

7.07 

7.39 

8.33 

9-99 

10.82 

10.94 

3.06 

9.88 

10.34 

9.86 

8.63 

8.00 

7.91 

12,500 

6.88 

7.10 

7.54 

8.75 

10.37 

10.90 

10.98 

3.17 

10.66 

10.95 

10.23 

9.03 

8.62 

8.57 

13,000 

6.86 

7.H 

7-71 

9.20 

10.62 

IO.96 

10.99 

8.70 

11.46 

11.55 

10.59 

9.54 

9-29 

9.26 

13,500 

6.91 

7.17 

7-95 

9-63 

10.77 

10.99 

11.01 

11.18 

12.25 

12.13 

11.00 

10.15 

10.00 

9-98 

14,000 

6.90 

7.23 

8.24 

10.01 

10.87 

11.00 

11.02 

12.11 

13.02 

12.69 

11.47 

10.83 

10.74 

10.73 

14,500 

6.90 

7-37 

8.66 

10.37 

10.94 

11.01 

11.02 

13.13 

13.79 

13-18 

11.98 

11.56 

11.51 

11.51 

15,000 

6.88 

7.44 

8.88 

10.52 

11.02 

11.02 

11.03 

13.97 

14.55 

13-81 

12.68 

12.32 

12.32 

12.31 


VJl 
U 0 
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TABLE VII.- PARTIAL COEFFICIENTS - Concluded 


ui 

■f=- 


T, 

°K 

(c) Prandtl number, Pr* 

(d) Levis number, Le' 

Pressure, atmospheres 

Pre s sure , atmo sphere s 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

100 

10 

1.0 

0.1 

0.01 

0.001 

0.0001 

2,000 

2,500 

3.000 

3.500 

4.000 

4.500 

5.000 

5.500 

6.000 

6.500 

7.000 

7.500 

8.000 

8.500 

9.000 

9.500 
10,000 

10.500 
11,000 

11.500 
12,000 

12.500 

0.775 

• 777 
.776 

• 771 
.764 
.753 
.744 
.740 

• 728 
.720 

• 711 
.698 
.691 
.687 
.687 
.695 
.704 
.710 

.716 

.719 

.719 

0.775 

.777 

.772 

• 758 

• 7^7 

• 739 
.732 
.724 
.715 
.700 
.690 
.690 
.693 
.707 
.714 
.725 

0.774 

• 773 

• 756 
.743 

• 739 

• 732 
.724 
.709 
.693 
.684 
.696 
.710 
.719 
.734 

0.773 

• 773 
.762 
.744 

• 738 
.732 
.723 
.705 
.687 
.689 
.703 
.710 
.721 

0.772 

.770 

.747 

.738 

.726 

.725 

.704 

.684 

.691 

.703 
• 711 

0.772 

.763 

.739 

• 734 

• 729 
.708 
.682 
.686 
.697 

0.772 

.750 

• 738 

• 734 
.717 
.683 
.683 
.694 

1.408 

1.445 

1.464 

1.449 

1.404 

1.327 

1.268 

1.235 

1.158 

1.118 

1.059 

• 979 
.893 
.805 
.726 
.662 
.618 
.585 
.564 

• 531 
.525 

1.407 

1.437 

1.424 

1.347 

1.265 

1.215 

1.167 

1.131 

1.071 

.975 

.860 

.756 

.669 

.618 

.584 

.558 

1.404 

1.405 
1.320 
1.225 
1.196 
1.158 
1.115 
1.019 

.885 
.746 
.652 
.601 
• 570 
.567 

1.355 

1.391 

1.327 

1.227 

1.180 

1.160 

1.105 

.982 

.808 

.669 

.605 

• 587 

• 570 

1.355 

1.355 

1.213 

1.167 

1.144 
1.114 
• 970 
.762 
.628 
.588 
.574 

1.346 
1.233 
1.157 
1.154 
1.130 
.998 
.750 
.610 
• 578 

1.335 

1.172 

1.147 

1.143 

1.051 

.786 

.609 

.580 


NACA TN 4150 


TABLE VIII.- CONVERSION FACTORS 


Multiply 

By 

To get 

In 

ZE ZH 

RT ° r RT 

0.0686 T (°K) 
0.287 T (°K) 
0.0686 T (°R) 
17l£ T (°R) 

Energy or enthalpy 

cal/gm 
joule /gm 
Btu/lb mass 
ft-lh/slug 

ZS ZCV 

TP R ' 

ZC ZC « 

R ’ R 

0.0686 

0.287 

0.0686 

1716 

Entropy or specific heat 

cal/gm °K 
joule/gm K 
Btu/lb mass R 
ft-lb/slug °R 

nAio 

1.462X10' 5 T 1/2 ^1 + 1 

(T in °K) 

2. 28x10" 8 T 1/2 (1 + 

(T in °R) 

Coefficient of 
viscosity 

gm/cm-sec 
lb sec /ft 2 

k k n 

ko or k Q 

4. 76 x 10" 6 t 1/2 ^1 + 
(T in °K) 

1.994X10" 5 t 1/2 ^1 + 
(T in °K) 
2.40x10'® T 1/2 ^1 + 
(T in °R) 

1.867X10' 3 T 1/2 (l + 
(T in °R) 

Coefficient of 
thermal conductivity 

cal 

cm- sec °K 

watt 
cm °K 

Btu 

ft-sec °R 

ft-lb 
ft-sec °R 


VJl 

VJl 
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Figure !•- Compressibility of air as a function of temperature. 
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Figure 2.- Energy of air as a function of temperature. 
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Figure 3 • - Entropy of air as a function of temperature . 
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Figure 4.- Specific heat of air at constant density as a function of 

temperature . 


Dimensionless speed of sound parameter, 
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Figure 5.- Zero frequency speed of sound parameter for air as a function of temperature. 
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Figure 7*- Stagnation temperature and pressure in air as a function of 

altitude and velocity. 


Coefficient of viscosity ratio, 


On 
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Figure 8 . - 


Ratio of the coefficient of viscosity for air to the reference coefficient, T} 0 , as a 

function of temperature. 
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Temperature, T, °K 


Figure 9 «- Ratio of the coefficient of thermal conductivity of air to the 
reference coefficient, k Q , as a function of temperature. 


Prandtl number, Pr 



Temperature, T, °K 

Figure 10.- Prandtl number for air as a function of temperature. 
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(a) Prandtl number 



(b) Lewis number 

Figure 11.- Partial coefficients for air as a function of temperature. 
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